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SECTION  I 


INTRODUCTION 


A.  BACKGROUND 

The  current  high  interest  in  the  effect  of  wind  shears  on  aircraft  has 
prompted  several  studies  (Refs.  1,  2,  and  5)  in  an  attempt  to  discover  and 
quantify  these  effects . 

Wind  shear  encounters  can  produce  catastrophic  results  by  introduction 
of  conditions  which  exceed  the  capability  of  the  pilot/vehicle  combination. 
Recent  examples  have  included  large  commercial  jet  transports  in  both 
takeoff  and  landing  flight  phases. 

The  purpose  of  this  study  was  to  explore  the  hazard  imposed  by  wind  shear 
on  different  classes  of  passenger  aircraft  and  to  consider  the  effects  of 
varying  operational  techniques  and  design  parameters . 

Two  techniques  were  used  to  study  the  problem.  First,  an  unmanned 
simulation  was  conducted  to  obtain  flight  profiles  in  wind  shear  conditions . 
Second,  an  analysis  was  performed  of  the  key  features  of  the  aircraft,  the 
pilot's  actions,  and  the  shears  themselves.  This  second  step  involved 
appropriate  linearization  of  pilot/vehicle  mathematical  models. 

The  unmanned  simulation  was  an  important  aspect  if  the  program.  The 
mathematical  model  which  took  the  role  of  the  real  pilot  was,  in  some 
respects,  superior.  If  we  assume  the  mathematical  model  to  be  valid,  it 
eliminated  pilot  variability  in  run- to-run  comparisons,  it  prevented 
learning  of  specific  shears,  and  it  was  less  expensive  than  using  real 
pilots . 

The  crucial  question  is,  of  course,  how  valid  is  the  mathematical  model 
pilot?  This  is  discussed  in  Section  II,  where  the  mathematical  pilot  model 
is  justified  as  being  representative  of  a large  class  of  human  pilots.  The 
pilot  model  was  not  designed  to  duplicate  exactly  a specific  pilot,  but 
instead,  to  approximate  "typical"  human  pilot  responses . To  help  to  insure 
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that  this  was  the  case,  the  pilot  model  used  in  this  study  was  checked 
against  and  "fine-tuned"  with  data  involving  human  pilots.  Additionally, 
studies  were  made  of  the  effect  of  variations  in  the  pilot  model  parameters 
to  check  for  any  overly  sensitive  aspects. 

This  study  specifically  involved  a comparison  of 

• Four  specific  aircraft 

(i)  A twin  engine  light  jet  transport  (referred  to 
as  the  "light  transport") 

(ii)  A wide-bodied  jumbo  jet  transport  (referred  to 
as  the  "jumbo") 

(iii)  A light  twin  engine  jet  (referred  to  as  the 
"business  jet"),  and 

(iv)  A STOL  commuter  aircraft  (referred  to  as 
the  "STOL  commuter") 

• Three  specific  wind  shears 

(ii  A nighttime  stable  boundary  layer  shear 
(ii)  A frontal  passage  shear 
(iii)  A thunderstorm  cold  air  outflow  shear 

• Both  approach  and  takeoff  cases 

• Variation  in  key  aircraft  performance  and  controllability 
parameters  and  all  pilot  model  parameters . 

Only  the  longitudinal  degrees  of  freedom  were  simulated  since  it  was 
believed  that  these  were  most  critical  in  determining  shear  hazard  and 
aircraft  safety. 

B.  REPORT  ORGANIZATION 

The  report  is  organized  into  the  following  five  sections  and  3 appendices . 
Section  II  contains  a description  of  the  pilot  model,  the  aircraft  models, 
the  shear  models,  and  the  simulation  procedure.  Section  III  contains  a 
description  of  the  linear  analysis  methods  and  procedures.  Section  IV 
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contains  the  results  from  the  approach  cases.  Section  V contains  the 
results  from  the  takeoff  cases.  Section  VI  presents  conclusions  and 
recommendations  based  on  this  simulation.  The  appendices  contain  back- 
ground information  concerning  the  pilot,  aircraft,  and  shear  models  used 
in  the  study. 
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SECTION  II 


SIMULATION  DESCRIPTION 


The  simulation  conducted  in  this  program  was  routine  except  for  one 
important  feature  — the  pilot  consisted  of  a mathematical  model.  In  the 
following  pages  the  rules  used  to  define  the  model  pilot's  control  strategy 
are  defined  and  discussed.  This  is  followed  by  descriptions  of  the  airplane 
model,  the  wind  shear  profiles,  and  the  simulation  procedure. 

A.  PILOT  MDDEL 

Three  distinct  forms  of  the  pilot  model  were  used.  These  were: 

(i)  A Conventional  Takeoff  and  Landing  (CTOL)  form 
for  use  with  the  light  transport,  jumbo,  and 
business  jet  aircraft  in  the  approach  flight 
phase, 

(ii)  A Short  Takeoff  and  Landing  (STOL)  form  for  use 
with  the  STOL  commuter  aircraft  in  the  approach 
flight  phase,  and 

(iii)  A takeoff  form  vised  for  all  aircraft  in  the 
takeoff  flight  phase. 

In  the  following  subsections  each  of  the  three  pilot  model  forms  will  be 
described  in  general;  however,  the  pitch  attitude  loop  which  is  common  to 
all  three  forms  will  be  described  only  in  connection  with  the  CTOL  pilot 
model.  Appendix  A presents  the  details  of  the  CTOL  pilot  model.  Each 
subsection  also  discusses  how  the  pilot  model  was  fine-tuned  to  match 
better  the  results  from  piloted  simulations. 

1 . The  CTOL  Pilot  Model 

An  overall  block  diagram  of  the  CTOL  pilot  model  is  given  in  Fig.  II- 1, 
the  key  features  of  which  include: 
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pe  Displacement  Error  Feedback 


Figure  II-1 . Block  Diagram  of  Pilot  Model  in  the  Context  of  Manual  Approach  Control 
Illustrating  Differences  in  Technique  for  CTOL  and  STOL 


• Regulation  of  pitch  attitude  with  elevator  command  via 
control  column, 

• Regulation  of  glide  slope  error  with  a commanded 
pitch  attitude, 

• Regulation  of  airspeed  error  with  throttle  with  a 5 kt 
threshold  of  indifference,  and 

• Crossfeed  from  pitch  command  to  throttle  for  large 
attitude  commands . 

These  four  key  features  are  based  on  the  quasi-linear  pilot  models  for 
compensatory  multi loop  systems  described  in  Ref.  4.  The  first  two  features 
have  been  extensively  analyzed  in  Refs . 5-8  in  the  absence  of  significant 
throttle  activity.  Validating  measurements  of  pilot  describing  functions 
representing  the  first  two  features  have  been  reported  in  Refs.  9— 11, 
again,  in  the  absence  of  significant  throttle  activity. 

The  basis  for  the  two  latter  features  are  not  nearly  as  well  documented 
or  validated  in  the  literature  on  the  subject  of  airspeed  regulation  and 
control  for  CTOL,  e.g.,  Ref.  5.  The  pilot  normally  inhibits  throttle 
activity  for  the  purpose  of  airspeed  regulation  with  a + 5 kt  threshold  of 
indifference  to  disturbance -induced  airspeed  variations  about  the  desired 
trimmed  airspeed.  If  the  aircraft  heave  damping-airspeed  product  (U  Z ) 
is  moderate  to  large  (e.g.,  as  in  the  cases  of  CTOL/jet  transports)  then 
only  small  transient  angle  of  attack  changes  need  be  provided  by  pitching 
control  in  order  to  regulate  glide  slope  deviation.  As  a result  of  the 
small  angle  of  attack  changes,  the  drag  fluctuates  only  slightly  and  the 
accompanying  airspeed  variations  are  usually  small.  In  contradistinction, 
if  U_Zw  is  small,  the  aircraft  can  experience  great  variations  in  speed. 

The  pilot,  however,  will  limit  his  commanded  variations  in  angle  of  attack 
(via  limits  on  commanded  pitch  attitude)  and  will  cross-feed  throttle 
commands  as  a function  of  his  desired  flight  path  angle  correction  in  order 
to  assist  his  otherwise  compensatory  airspeed  regulation  via  throttle 
(Ref.  5). 

The  downside  pitch  attitude  limit,  9^  , in  the  pilot  model  (Fig.  II-1 
and  Fig.  A-1 ) was  set  at  a value  which  would  give  the  pilot  model  a suffi- 
cient amount  of  nose  down  capability  but  would  not  allow  it  to  generate 
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unrealistically  large  negative  6 values.  Since  the  shear  profiles  used 
primarily  drive  the  aircraft  low  and  slow  and  result  in  a positive  A0, 

9min  WHS  no^  a parameter  and,  indeed,  was  never  exceeded  during 

the  runs . However,  the  upper  pitch  attitude  limit,  0 , was  felt  to  be 

an  important  parameter  since  hitting  this  limit  caused  the  pilot  model  to 
apply  and  hold  full  power,  i.e.,  to  essentially  abort  the  approach  and 
initiate  a go- around. 

It  should  be  noted  that  the  function  of  the  limit  was  to  prohibit 

the  pilot  model  from  placing  the  aircraft  in  a pitch  attitude  that  was 
unrealistically  large,  i.e.,  a pitch  attitude  that  a human  pilot  would  be 
reluctant  to  use.  The  reasoning  behind  having  the  pilot  model  apply  and 
hold  full  power  if  0 was  exceeded  was  that  the  pilot,  having  used  al 1 
of  his  available  pitch  authority,  would  then  use  the  only  other  control 
available,  thrust,  in  an  attempt  to  extricate  himself  from  the  situation. 

Several  methods  of  setting  0 were  considered.  The  method  finally 

max  “ 

used  was  to  set  0 equal  to  the  steady  state  pitch  attitude  which  would 
max 

be  achieved  at  the  designated  approach  speed  with  maximum  thrust. 

The  CTOL  pilot  model  was  fine-tuned  by  comparing  its  responses  for  the 
light  transport  with  two  manned  simulator  approaches  until  the  responses 
of  the  pilot  model  compared  favorably  prior  to  and  throughout  the  recovery 
(past  the  peak  airspeed  and  altitude  deviations)  from  the  wind  shear  upset. 
The  timing  and  magnitudes  of  the  responses  of  the  aircraft  controlled  by 
the  pilot  model  were  acceptably  similar  to  those  from  the  manned  cases. 

2.  The  STOL  Pilot  Model 

The  STOL  pilot  model  (Fig.  II-2)  is  conceptually  similar  to  the  CTOL 
version  (Fig.  II-1 ) except  that  the  STOL  version  incorporates 

• A glide  slope  error  to  throttle  feedback  without  a 
threshold 

• An  airspeed  error  to  pitch  attitude  feedback  with  a 
5 kt  threshold 

• No  crossfeeds. 
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STOL  pilot  models  have  been  extensively  analyzed  in  Refs.  12-15* 

Results  of  several  manned  simulations  designed  to  validate  STOL  piloting 
techniques  are  reported  in  Refs.  15-20.  A number  of  STOL  flight  test 
validations  are  reported  in  the  bibliography  of  Ref.  20,  which  provides  a 
concise  summary  and  comparison  of  CTOL  and  STOL  piloting  techniques . 

The  attitude  compensation,  Y0,  for  the  STOL  pilot  model  was  determined 
in  the  same  manner  as  for  the  CTOL  pilot  model.  The  variation  of  glide 
slope  loop  gain  with  range  is  identical  to  that  in  the  CTOL  pilot  model. 

A detailed  list  of  the  STOL  pilot  model  parameters  is  presented  in  Appendix  A, 
Table  A-1 . 

As  with  the  CTOL  model,  the  STOL  pilot  model  was  fine-tuned  by  comparison 
with  piloted  simulation  results.  In  this  case,  however,  the  piloted  simu- 
lation results  were  extremely  limited,  and  caution  should  be  exercised  in 
applying  this  pilot  model  to  any  STOL  aircraft  other  than  that  for  which  it 
was  adapted.  The  piloted  data  that  were  available,  however,  were  matched 
utilizing  the  same  phase  margins  and  crossover  frequencies  for  the  glide 
slope  and  airspeed  loops  as  in  the  CTOL  pilot  model.  No  requirement  for 
crossfeeds  was  apparent  from  the  comparison  with  the  piloted  simulation 
results . 

3.  Takeoff  Pilot  Model 

The  takeoff  pilot  model  involves  only  airspeed  and  pitch  attitude  feed- 
back loops  in  series  to  the  elevator  (Fig.  II-3).  This  model  was  used  for 
all  aircraft.  For  takeoff,  thrust  is  left  unchanged  at  its  initial  takeoff 
setting,  and  speed  variations  are  controlled  with  pitch  attitude. 

It  was  assumed,  in  the  absence  of  other  information,  that  the  same 
pitch  attitude  and  airspeed  loop  phase  margins  and  crossover  frequencies 
would  apply  for  the  takeoff  case  as  for  the  approach  cases.  The  only 
alteration  which  was  made  was  the  speed  threshold,  which  was  reduced 
from  5 kt  to  3 kt. 

It  is  worthwhile  to  note  that  for  the  takeoff  case  the  effect  of  the 

shears  is  to  retard  the  speed  of  the  aircraft,  thus  resulting  in  a nose-down 

pitch  input.  Thus,  0 . becomes  more  critical  than  0 in  this  pilot 
’ min  max 
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Figure  II-?.  Takeoff  Pilot  Model  (All  Aircraft) 


model.  No  rational  justification  for  choosing  a other  than  pilot 

reluctance  to  lower  the  nose  further  was  suggested,  and  a 0^  of  -10  deg 
was  chosen  arbitrarily. 

A detailed  list  of  the  takeoff  pilot  model  parameters  for  the  cases 
and  aircraft  tested  is  presented  in  Appendix  A,  Table  A-2. 

B.  AIRCRAFT  MODEM 

To  facilitate  the  modeling  of  a variety  of  different  aircraft  for  this 
simulation,  a simple  generic  aircraft  model  form  was  devised.  The  simplicity 
of  this  model  allowed  rapid  modeling  of  each  aircraft  and  decreased  the  time 
required  to  check  each  aircraft,  yet  retained  the  important  aerodynamic, 
control  system,  and  propulsion  characteristics  of  each  individual  aircraft. 

Four  aircraft,  each  representing  a distinct  class  of  commercial  trans- 
ports, were  simulated.  The  aircraft  are  described  briefly  below: 

• The  light  transport  — A twin-engine,  low- wing,  light  , 
jet  transport  with  a takeoff  gross  weight  of  110,000  lb. 

* 

• The  jumbo  — A four- engine,  low- wing,  jumbo  jet  with  * 
a takeoff  gross  weight  of  753,000  lb. 

k-» 

• The  business  jet  — A twin-engine,  light,  swept-wing, 
business  jet  with  a takeoff  gross  weight  of  12,500  lb. 

• The  ST0L  commuter  — A twin-engine,  turbo  prop, 
commuter  airline  ST0L  aircraft  with  a takeoff  gross 
weight  of  1 1 , 579  lb . 

Each  model  has  been  used  previously  in  manned  simulator  experiments  at 
NASA  Ames  Research  Center  (ARC ) . > 

% 

The  jumbo  and  light  transport  models  were  developed  by  their  manufacturer 
based  on  the  data  from  Refs.  21  and  22.  The  business  jet  model  was  originally 
developed  by  STI  for  a vortex  encounter  study  and  is  described  in  Ref.  23. 
Finally,  the  ST0L  commuter  model  originated  at  NASA  ARC  for  guidance  and 
navigation  research  and  is  described  in  Ref.  2k. 

Each  of  the  above  models  utilizes  a table  look-up  scheme  for  basic 
aerodynamic  and  propulsion  features.  Lift  and  pitching  moment  are  described 
to  a point  beyond  stall,  and  engine  thrust  response  varies  appropriately 
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with  thrust  level.  Control  system  descriptions  include  servo  lags  and 
limits  on  surface  rate  and  displacement. 


Basic  aircraft  design  features  are  summarized  in  Table  II-1  and  in  the 
general  arrangement  drawings  in  Fig.  Xl-t.  Stability  derivatives  and  trans- 
fer functions  are  listed  in  Appendix  B using  the  notation  defined  in 
Ref.  25. 

C.  SHEAR  MODELS 

Three  different  wind  shear  models  were  provided  by  NASA  for  use  in  this 
study.  They  are  identified  as: 

• Boundary  layer  shear 

• Logan  shear 

• Kennedy  shear . 

Each  is  defined  in  terms  of  a "look-up"*  table  of  wind  velocity  (horizontal 
and  vertical  components)  versus  altitude,  although  a range  dependency  was 
also  contrived  for  the  latter  two  shears.  (The  dependency  aspect  is  dis- 
cussed fully  at  the  beginning  of  Section  IV.)  The  three  shears  are  summarized 
in  Fig.  II-5. 

One  shear  represents  a nighttime  stable  boundary  layer  shear  corresponding 
to  a low  level  temperature  inversion  overlayed  by  fairly  strong  winds 
immediately  above  the  inversion.  The  simulated  velocity  profile  is  given 
in  Table  II-2. 

The  Logan  shear  represents  a shear  caused  by  a frontal  passage.  It 
consists  of  a 33  kt  tailwind  at  1600  ft  which  changes  gradually  to  a 14  kt 
headwind  at  200  ft  then  rapidly  to  a L kt  tailwind  at  20  ft.  The  altitude 
and  velocity  breakpoints  are  given  in  Table  II-3-  The  critical  aspect  of 
this  shear  is  the  decreasing  headwind  below  200  ft.  Nominally,  the  altitude 
dependent  version  of  this  shear  should  cause  greater  aircraft  problems  than 


* A "look-up"  table  implies  the  use  of  an  interpolation  algorithm  in  the 
computer  simulation. 
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TABLE  II- 1 

SUMMARY  OF  AIRCRAFT  MODEL  CHARACTERISTICS 


JUMBO 

LIGHT 

TRANSPORT 

BUSINESS 

JET 

STOL 

COMMUTER 

Wing  Arsa  (ft2) 

5,500 

980 

231 .77 

420 

Wing  Aspect  Ratio 

6.96 

8.83 

5-46 

10.06 

Landing 

Weight  (lb) 

570,000 

95,000 

12,500 

11,000 

Flap  Deflection  (deg) 

30 

4o 

40 

40 

Cltoax 

2.44 

3.0 

1 .61 

2.75 

Takeoff 

Weight  (lb) 

733,000 

110,000 

12,500 

11,579 

Flap  Deflection  (deg) 

10 

5 

20 

20 

1 .86 

2.0 

1.48 

2-53 
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TABUS  II-2 


BOUNDARY  IAYER  SHEAR 


HORIZONTAL  WIND 

ALTITUDE,  h (:Dt)  VEIflCITY,  ug*  (kt) 


1500 

-22 

1300 

-22 

1200 

-23 

1050 

-23 

900 

-25 

750 

-26 

600 

-30 

450 

-32 

300 

-3- 

C\J 

1 

150 

-20 

75 

-12 

20 

-8 

A positive  u corresponds  to  a tailwind. 

o 

The  vertical  gust  component,  w so. 

6 


TABLE  II-3 


LOGAN  SHEAR 


ALTITUDE,  h (ft) 

1600 

1100 

900 

800 

700 

600 

500 

400 

300 

200 

120 

20 

0 


HORIZONTAL  WIND 
VELOCITY,  u * (kt) 

■■  ■ Si 

53 

26 

23 

22 

21 

20 

17 

2 

-11 

-14 

-5 

4 

4 


# 


A positive  u corresponds  to  a tailwind. 
6 

The  vertical  gust  component,  w ■ 0 . 

S 
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the  range -dependent  version,  since  the  shear  causes  the  aircraft  to  go  low 
and  slow  which,  for  an  altitude -dependent  shear,  causes  further  decrease 
in  headwind. 

The  Kennedy  shear  is  the  most  complex  of  the  shears  with  both  horizontal 
and  vertical  components . The  shear  is  representative  of  the  conditions 
due  to  a thunderstorm  cold  air  outflow  and  was  derived  from  an  analysis  of 
the  conditions  encountered  by  Eastern  Airlines  Flight  066  near  John  F. 

Kennedy  Airport  in  June  1 975-  The  wind  profiles  are  shown  in  Table  Il-k. 
Initially  the  aircraft  encounters  a mild  headwind.  As  the  aircraft  descends, 
it  encounters  a simultaneous  increase  in  headwind  and  a downdraft  (near 
800  ft).  Near  I50  ft  the  headwind  starts  to  decrease  and  the  downdraft 
changes  briefly  to  an  updraft.  The  winds  become  particularly  adverse 
around  300  ft  where  there  is  a downdraft  and  a very  rapidly  decreasing 
headwind. 

The  Logan  and  Kennedy  shears  were  investigated  from  both  range-dependent 
and  altitude-dependent  points  of  view  for  reasons  we  shall  discuss  at  the 
beginning  of  Section  IV.  For  the  altitude- dependent  shears,  the  wind  velocity 
look-up  table  was  entered  with  the  vehicle  altitude,  h^,  while  for  the  range- 
dependent  shears,  the  look-up  table  was  entered  with  the  altitude  of  the 
glide  slope,  h_e,  corresponding  to  the  vehicle's  ground  range  from  the  glide 
path  intercept  point  (GPIP)  sketched  in  Fig.  II-6.  The  differences  between 
the  range-  and  altitude-dependent  versions  of  a shear  are  small  if  the 
aircraft  stays  close  to  the  glide  slope,  but  can  be  significant  if  large 
deviations  from  the  glide  slope  occur.  The  mechanization  of  the  shears 
indicated  above  was  used  for  both  the  approach  and  takeoff  flight  phases. 

For  the  takeoff  the  nominal  flight  path  angle  used  to  compute  the  wind 
velocities  for  range- dependent  shears  was  3 deg.  For  the  altitude- dependent 
shears  the  height  of  the  aircraft's  eg  was  used. 

Since  the  aircraft  normally  had  a takeoff  climb  gradient  greater  than 
3 leg,  the  altitude-dependent  cases  usually  caused  more  rapid  changes  in 
the  wind  velocity  than  did  the  range -dependent  cases . However,  if  the 
shears  caused  the  climb  gradient  to  be  less  than  3 deg,  the  altitude- 
lependent  cases  would  cause  less  rapid  wind  changes  than  the  range-dependent 
cases . 
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TABLE  II-4 


kennt;d>  n hear 


Horizontal  Component 


WIND  VELOCITY, 


ALTITUDE,  h (ft) 

u * (kt) 

6 

1500 

8 

750 

8 

450 

30 

300 

15 

200 

-4 

0 

-10 

Vertical  Component 

WIND  VELOCITY, 

ALTITUDE,  h (ft) 

wg*  (kt) 

850 

0 

800 

7 

700 

9 

600 

10 

500 

8 

400 

-10 

300 

4 

250 

3.5 

200 

3 

0 

0 

* 


A positive  u corresponds  to  a tailwind  and 

o 

a positive  w , a downdraft. 

6 
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D.  SIMULATION  PROCEDURE 


The  standard  procedure  used  to  initiate  a simulated  approach  run  was 
to  position  the  aircraft  at  an  altitude  of  1(500  ft  on  speed  and  on  glide 
slope.  For  the  CTOL  aircraft  a 3 deg  glide  slope  was  used.  The  STOL 
commuter  aircraft  was  nominally  flown  on  a 6 deg  glide  slope,  but  some 
runs  also  were  made  with  a 3 deg  glide  slope.  The  aircraft,  initially 
trimmed,  was  then  allowed  to  fly,  under  the  control  of  the  pilot  model, 
down  the  glide  slope  and  through  the  shear.  The  run  was  automatically 
terminated  whenever  the  aircraft  hit  the  ground  or  overflew  the  glide  slope 
transmitter.  The  values  for  the  nominal  approach  speeds,  sink  rates,  flight 
path  angles,  and  landing  weights  for  the  aircraft  are  given  in  Table  II-5. 

For  the  takeoff  tests,  the  aircraft  heading  relative  to  the  horizontal 
wind  was  reversed.  The  aircraft  would  be  positioned  at  an  altitude  of 
30  ft  at  a longitudinal  position  that  would  place  the  aircraft  outbound  on 
a 3 deg  glide  slope,  regardless  of  the  actual  takeoff  climb  angle  of  the 
aircraft  (see  Fig.  II-7).  This  was  done  so  that,  regardless  of  the  initial 
climb  angle  of  the  aircraft  being  investigated,  the  wind  conditions  at  the 
30  ft  trim  point  would  be  identical  for  the  range-  and  altitude- dependent 
versions  of  a given  shear. 

The  aircraft  was  trimmed  at  this  initial  position  in  the  presence  of 
the  appropriate  winds  at  a particular  speed  and  power  setting.  It  then 
proceeded  to  fly,  under  the  guidance  of  the  pilot  model,  through  the  shear. 
The  nominal  takeoff  speeds,  climb  rates,  flight  path  angles,  and  takeoff 
weights  for  each  aircraft  are  listed  in  Table  II-6. 

As  the  aircraft  proceeded  along  its  flight  tisrough  the  shear,  the 
following  variables  were  recorded  on  strip  charts: 

• Aircraft  eg  height  above  ground 

• Glide  slope  error 

• Sink  rate 

• Vertical  deviation  from  glide  slope 

• Horizontal  wind  component 
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Glide  Slope 


Figure  II-6.  Entry  Altitudes  for  Wind  Velocity  Tables 
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TABLE  II-5 

NOMINAL  APPROACH  CONDITIONS 


AIRCRAFT 

APPROACH 

SPEED 

(kt) 

Jumbo 

142 

Light  Transport 

130 

Business  Jet 

125 

STOL  Commuter 

70 

SINK  RATE 
(ft/sec) 

FLIGHT 

PATH 

ANGLE 

(deg) 

landing 

WEIGHT 

“(lb) 

-12.6 

-3 

570,000 

-11  .5 

-3 

95,000 

-11  .0 

-3 

12,500 

-12.4 

-6 

11,000 
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I 


[ 


7ac  = aircraft  climb  angle 


Figure  II-7 • Takeoff  Initial  Conditions 
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TABLE  II-6 


NOMINAL  TAKEOFF  CONDITIONS 


AIRCRAFT 

TAKEOFF 

SPEED 

CLIMB  RATE 

FLIGHT 

PATH 

ANGLE 

TAKEOFF 

WEIGHT 

(kt) 

(ft/sec) 

(deg) 

(Id) 

J umbo 

187 

38.5 

7.0 

733,000 

Light  Transport 

161 

37-^ 

7-9 

110,000 

Business  Jet 

135 

62.5 

15-9 

12,500 

STOL  Commuter 

78 

21  .3 

9-3 

1 1 >579 

• Throttle  position 

• Engine  thrust,  lb 

• Pitch  attitude 

• Vertical  wind  component 

• Engine  thrust,  $ max 

• Elevator  position 

• Angle  of  attack 

• Velocity  error 

• Aircraft  range  from  glide  slope  transmitter 

• Various  internal  pilot  model  signals  and  flags 
indicating  the  internal  pilot  status. 

In  addition  to  the  above  strip  chart  variables,  the  digital  computer 
was  programmed  to  monitor  eight  variables  and  record  the  maximum  and  minimum 
value  of  each  variable,  along  with  the  time  and  altitude  at  which  the 
maximum  and  minimum  values  occurred.  These  eight  variables  were: 

• Rate  of  climb 

• Pitch  angle 

• Angle  of  attack 

• Engine  pressure  ratio 

• Column  position 

• Glide  slope  error  (angular) 

• Vertical  distance  from  glide  slope 

• Velocity. 

Examples  of  typical  digital  output  are  given  in  Table  II-7. 
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SECTION  III 


LINEAR  ANALYSIS 


A.  BACKGROUND 

The  pilot  and  aircraft  models  used  in  the  simulation  effort  included 
nonlinearities.  The  pilot  model  involved  thresholds  and  limiters  while 
the  aircraft  models  incorporated  nonlinear  aerodynamic,  control  system, 
and  propulsion  system  effects . These  nonlinearities  may  be  required  in 
a simulation  to  achieve  a valid  and  reasonable  correlation  with  real  world 
results,  but  they  tend  to  obscure  key  effects  and  parameters.  To  counter 
this  problem,  linearized  models  were  developed  to  simplify  the  pilot  and 
vehicle  and  reveal  the  essential  features  of  both. 

The  use  of  linear  small  perturbation  aircraft  models  allows  the  expres- 
sion of  vehicle  dynamic  characteristics  with  a minimum  number  of  essential 
parameters.  Additionally,  the  linear  models  contain  explicitly  the  vehicle 
stability  derivatives  which  have  a key  effect  on  aircraft  behavior.  Ex- 
pressed in  this  way,  the  key  dynamic  parameters  of  the  model  are  more 
clearly  partitioned.  The  stability  derivatives  have  the  additional  benefit 
that  each  one  is  easily  derived  from,  or  associated  with,  an  aerodynamic  or 
geometric  characteristic  of  the  aircraft. 

In  the  following  pages  we  shall  define  the  linear  models  used  in  our 
analysis  and  proceed  to  lose  these  models  in  exposing  dominant  wind  shear 
effects . We  shall  end  this  section  by  presenting  a simplified  view  of  the 
mechanism  of  wind  shear  acting  on  the  piloted  aircraft.  This  will  serve 
as  a background  for  the  subsequent  presentation  of  simulation  results. 

B.  LINEARIZED  PILOT  AND  AIRCRAFT  MODELS 

The  pilot  models  were  linearized  by  removing  the  thresholds  and  limits 
from  all  pilot  models  and  by  removing  the  crossfeed  term  from  the  CTOL 
approach  pilot  model.  Figures  III-1 , III-2,  and  III-3  show  block  diagrams  of 
the  linear  pilot  models  for  the  CTOL  approach,  STOL  approach,  and  takeoff, 
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Figure  III-1 . Linearized  CTOL  Approach  Pilot  Model 


Figure  III-2.  Linearized  STOL  Approach  Pilot  Model 


Figure  III-3.  Linearized  Takeoff  Pilot  Model 


respectively.  The  numerical  values  of  transfer  functions  correspond  to 
those  of  the  non-linear  model  (see  Appendix  A). 

The  aircraft  models  were  linearized  by  the  conventional  methods  of 
small  perturbations  about  a nominal  trim  point  and  expressed  in  terms  of 
dimensional  stability  derivatives  in  a body-fixed  stability  axis  system. 
Both  the  aircraft  aerodynamic  and  propulsion  system  characteristics  were 
treated  in  this  matter.  Reference  25  may  be  cons  lilted  for  details  of  the 
equations  of  motion.  The  stability  derivatives  for  the  models  at  base- 
line conditions  are  given  in  Appendix  B. 

The  linearized  descriptions  of  pilot  and  aircraft  dynamics  can  be 
directly  utilized  by  formulating  closed  loop  transfer  functions  of  air- 
craft motion  due  to  applied  gust  disturbances.  As  we  shall  discuss  in 
Section  IV,  one  of  the  most  significant  motion  quantities  is  altitude 
excursion.  Hence,  it  is  instructive  to  consider  the  transfer  functions: 


Ah  , Ah 

— and  — 

U„  Wp 


Two  examples  of  the  above  transfer  functions  are  shown  in  Figs.  Ill-t- 
and  III-5  for  the  light  transport  and  in  Figs.  III-6  and  III-7  for  the  STOL 
commuter.  The  light  transport  is  representative  of  all  the  CTOL  types 
considered  in  this  study.  Some  important  general  observations  can  be  made 
from  these  figures.  First,  the  frequency  response  characteristics  for  both 
types  of  aircraft  are  comparable,  although  the  STOL  response  is  not  as  well 
damped  as  the  CTOL  response.  The  pilot-vehicle  combination  acts  as  a band 
pass  filter  to  horizontal  gust  velocity  and  as  a low  pass  filter  to  vertical 
gust  velocity.  (We  can  also  characterize  the  horizontal  gust  response  as  a 
low  pass  filter  for  horizontal  gust  acceleration  or  rate  of  shear,  i.e., 
u .)  One  of  the  main  features  of  the  figures  is  that  closed- loop  pitch 
response  modes  which  are  characteristically  above  1 rad/sec  are  well 
attenuated.  The  predominant  aircraft  response  occurs  in  the  spectral 
range  below  1 rad/sec.  Further,  by  neglecting  the  pitch  response  and 
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.!□  “ (rad/sec) 


Closed-Loop  Altitude  Response  to  Horizontal  Gust  Velocity  Component 
(STOL  Commuter,  Approach) 


r 


concentrating  on  the  predominant  path  response,  we  can  simplify  the  dynamics 
enough  to  identify  easily  the  key  mechanisms  of  wind  shear  effects  on  air- 
craft motion. 

C.  SIMPLIFIED  PILOT- VEHICUS  DYHAMIC8 

The  key  to  describing  the  predominant  effects  of  wind  shear  on  aircraft 
motion  is  to  utilize  pitch- cons trained  equations  of  motion.  In  effect,  we 
assume  that  the  pilot  can  instantaneously  obtain  any  pitch  attitude  he 
commands,  and  that  gusts  have  no  effect  on  pitch  attitude  (i.e.,  pitch 
attitude  is,  itself,  an  independent  control).  These  are  not  valid  assump- 
tions if  one  is  concerned  about  pitching  motion,  per  se,  but  they  are 
reasonably  valid  if  flight  path  and  airspeed  are  the  main  concern  as  in 
the  case  here.  Further,  pitching  motion  can  be  handled  as  a second  order 
effect  within  the  simplified  context,  if  desired. 

The  detailed  derivation  of  simplified  pilot-vehicle  dynamics  is  carried 
out  in  Appendix  C.  The  starting  point,  however,  is  worth  stating  here  in 
terms  of  the  following  equations  of  motion: 


S - xu 

Xws 

ua 

• 

(*a  “ 

6)  X§rji  — S 

-Xw 

Zu 

s(s  - Z„) 

Ah 

;za 

— Z&iji  o 

Zw 

with 


&T  = — KuUa 
= -K^ 


0 = 


Pilot's  CTOL 
technique 


0 

&T 

ug 


w, 


8 


or 


Sf  = — i 


0 = — ] 


K(j^h 

'^uua 


Pilot's  STOL 
technique 


The  first  benefit  of  this  description  is  that  we  have  vastly  reduced  the 
number  of  essential  variables  which  describe  even  the  linearized  vehicle. 
Note  that  the  airspeed  and  flight  path  motion  due  to  attitude,  throttle, 
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horizontal  gusts,  and  vertical  gusts  is  dependent  only  upon  the  six  parameters 
Xu>  Xw,  Zjj,  Z w,  V,  and  = V ^ and  = V Z^).  Each  of  these  is 

easily  estimated  as  demonstrated  in  Ref.  20.  Two  pilot  parameters  Eire  in- 
volved, Ku  and  K4.  These  are  set  directly  in  proportion  to  their  respective 
control  loop  tightness  as  expressed  by  crossover  frequencies,  co^  and  o^. 

For  the  CTOL  case; 


and 


~ Kd^u  = caC(i 


The  altitude  response  due  to  horizontal  and  vertical  gusts  can  be 
expressed  in  the  following  generic  forms: 


or 


and 


Ah 

ug 

Ah 

Ah 

Wg 


_ 2g 


(s  + J-)(s2  + 2?^ds  + cn|) 


Lu 


_ 2S 
V 
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, 1 X*  . 
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w T u Z u7 
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where 


and 


Tu 


2 

“d 
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Zw  represents  heave  damping  and  is  equal  to 
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These  expressions  can  be  sketched  in  terms  of  amplitude  response  to  horizontal 

gust  velocity,  u , horizontal  gust  acceleration  or  rate  of  shear,  ut,  and 
6 

vertical  gust  velocity,  w , as  shown  in  Fig.  III-8.  Accordingly,  altitude 

6 

response  to  horizontal  wind  shear  is  inversely  proportional  to: 

• Airspeed 

• Airspeed  loop  tightness  (<uCu) 

• Square  of  flight  path  loop  tightness  (o>C(i) 

for  the  spectral  region  extending  out  to  the  flight  path  regulation  mode 
at  frequency  Altitude  response  to  vertical  gust  velocity  is  proportional 

to: 


• Airspeed 

• Lift  curve  slope, 

and  inversely  proportional  to: 

• Wing  loading 

• Square  of  flight  path  loop  tightness  (o>„.). 

These  results  will  be  useful  for  aiding  in  the  interpretation  of  simula- 
tion results  to  be  presented  in  the  next  section. 


ensitivity  to  u (gust  rate) 


log 


Sensitivity  to  w 


& 


log 


Figure  III-8.  Sketches  of  Altitude  Amplitude  Response  Asymptotes 
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SECTION  IV 


APPROACH  RESULTS 


A.  BACKGROUND 

The  presentation  of  wind  shear  hazard  results  for  the  approach  flight 
phase  is  divided  into  three  parts  with: 


• Baseline  aircraft  and  pilot  models 

• Variation  of  operational  techniques  which  affect 
shear  hazard 

• Variation  of  propulsion  system  features  which 
affect  shear  hazard. 


Prior  to  discussing  these  three  parts  of  the  results,  however,  it  is  im- 
portant to  define  the  metric  used  to  describe  wind  shear  hazard  and  certain 
aspects  of  modeling  wind  shear  which  impact  the  interpretation  of  the 
results . 


Several  metrics  were  considered  including  airspeed  loss  and  minimum 
altitude  reached  prior  to  recovery  from  the  disturbance.  The  most  direct 
and,  we  believe,  the  most  meaningful  proved  to  be  maximum  vertical  excursion 
below  the  glide  slope  or  nominal  flight  path.  This  is  illustrated  in 
Fig.  IV- 1 . 

Initially,  airspeed  loss  was  considered  as  a candidate,  and  although 
airspeed  loss  was  frequently  substantial,  its  use  was  ultimately  discarded 
for  two  reasons.  First,  airspeed  loss  is  mainly  a function  of  shear  abrupt- 
ness. The  features  which  would  have  the  most  direct  impact  on  airspeed  loss 
are  the  convective  (or  speed  damping)  tendency  of  the  bare  airframe  and  the 
speed  regulation  abilities  of  the  pilot.  This  relationship  can  be  easily- 
expressed  as: 


-u 


u 


s + 


u 
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where  -u 
a 


airspeed  loss 


u 

g 


A 

A 

A 


tailwind  increase 


time  constant  representing  the  aircraft's 
convective  tendency  and/or  the  pilot's 
speed  regulation  ability 


The  above  relationship  implies  that  airspeed  loss  matches  the  tailwind 
increase  for  a period  nearly  equal  to  T . At  best,  T^  includes  time  to 
recognize  an  airspeed  change,  to  take  effective  control  action,  and  for  the 
control  itself  to  produce  a longitudinal  acceleration.  This  could  typically 
be  5 to  10  seconds  based  on  observations  from  manned  simulations.  Never- 
theless, the  airspeed  loss,  per  se,  is  directly  tied  to  the  shear  itself. 

The  second  shortcoming  of  airspeed  loss  as  a metric  is  its  indirect 
relation  to  angle  of  attack.  Airspeed  loss  does  not  directly  result  in  a 
large  angle  of  attack  increase.  During  horizontal  shear  encounters,  the 
angle  of  attack  remains  relatively  constant  as  the  airspeed  decreases,  with 
only  3 to  5 deg  changes  in  a observed.  Also,  during  vertical  shear  encounters, 
the  sink  rate  increases  without  any  change  in  airspeed. 

While  we  recognize  the  value  of  airspeed  and  angle  of  attack  as  legiti- 
mate safety  margin  indicators,  we  claim  that  the  ultimate  metric  of  hazard 
should  be  a precursor  of  ground  contact  itself.  Hence,  we  are  forced  to 
consider  an  altitude-related  quantity. 

One  example  of  an  altitude  metric  is  the  absolute  minimum  altitude 

encountered  during  recovery  from  a shear  disturbance.  This  is  meaningful, 

however,  only  if  a local  minimum  occurs  in  altitude  (it  may  not  occur)  and 

if  the  excursion  is  tied  to  an  absolute  altitude  of  shear  onset.  Such  a 

set  of  conditions  is  unnecessarily  restrictive.  Thus,  we  return  to  the  simple 

yet  effective  Ah  metric  defined  in  Fig.  IV-1 . 

J max 

Another  matter  which  requires  discussion  prior  to  the  presentation  of 
results  is  the  way  the  wind  shear  is  modeled.  Despite  the  use  of  actual 
wind  shear  encounters  as  the  basis  of  this  study,  the  results  can  vary 
depending  upon  how  one  implements  the  shear  profiles . 
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An  actual  shear  profile  is  a function  of  spatial  position  (i.e., 
range  and  altitude)  and  time.  Wind  shear  measurements  based  upon  a single 
trajectory  through  an  air  mass  cannot  have  the  full  dimensionality  required, 
and  an  arbitrary  dependency  must  be  assigned.  It  is  important  to  recognize 
that  shear  profiles  with  differing  kinds  of  dependencies  will  produce  dif- 
ferent simulator  results.  This  is  illustrated  in  Fig.  IV-2. 

Let  us  further  consider  one  of  the  kinds  of  shear  dependencies. 

Figure  IV-3  shows  the  effect  of  an  altitude- dependent  horizontal  wind  shear 
on  an  attitude-stabilized  pilot-vehicle  combination.  If  the  wind  shear  is 
approximately  linear  with  altitude,  a linear  feedback  occurs . In  the  case 
of  a decreasing  headwind  the  effect  is  destabilizing.  In  fact  if  the  change 
in  headwind  with  altitude  is  approximately  greater  than  g/v  the  total  pilot- 
vehicle-shear  combination  is  unstable. 

Range -dependent  shears  create  a similar  but  far  weaker  closed  loop 
effect.  Time- dependent  shears  are  essentially  open  loop  (i.e.,  independent) 
disturbances . 

In  -view  of  the  foregoing,  we  therefore  studied  the  effects  of  wind  shear 
using  both  an  altitude- dependent  and  a range-dependent  model  of  the  Kennedy 
and  Logan  shears . However,  only  an  altitude- dependent  boundary  layer  shear 
was  used. 

It  will  be  apparent  that  the  altitude- dependent  shears  are  usually  more 
severe  than  their  range -dependent  counterpart.  In  view  of  the  insufficient 
knowledge  of  the  real  shear  dependency,  we  can  only  speculate  that  the  true 
degree  of  severity  is  somewhere  between  the  range-dependent  and  altitude- 
dependent  solutions. 

B.  RESULTS  FOR  BASELINE  AIRCRAFT  AID  PILOT  MDDELB 

Figure  IV- 4 depicts  the  maximum  deviation  below  the  glide  slope,  Ah^^, 
for  the  four  aircraft  during  encounters  with  the  boundary  layer  shear  and 
the  altitude- dependent  versions  of  the  Logan  and  Kennedy  shears.  (Note  that 
the  boundary  layer  shear  was  used  only  with  the  light  transport  and  STOL 
commuter.)  The  data  for  the  range -dependent  Logan  and  Kennedy  shears  are 
presented  in  Fig.  IV-5.  The  data  indicate  that,  for  all  aircraft,  the 
shears  can  be  rated  in  order  of  increasing  hazard: 
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TIME  DEPENDENT  SHEAR 


TIME 


Eigure  IV-f! 
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1 Open  loop  interaction  between 
shear  and  pilot-vehicle  "I 


RANGE  DEPENDENT  SHEAR 


Closed  loop  interaction  between 
shear  and  pilot- vehicle ) 


APTITUDE  DEPENDENT  SHEAR 
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Varieties  of  interactions  Between  Wind  Shear  and  Pilot-Vehicle 
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TKJTE:  TIE  PIL0T-VEHICLE-S1IEAR  COMBINATION  IS  UNSTABLE  FOR 

Sus<ILlS  • i 

Hi  , " v 

'u 

du 

e.g..  if  V = 120  kt.  the  critical = 9.4  ___  headwind- to- tailwind) 

figure  IV *o.  The  Effect  of  an  Altitude-Dependent  Horizontal  Wind  Shear 
on  the  Attitude-Stabilized  Pilot- Vehicle  Combination 
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BOUNDARY  LAYER  3 HEAR 


LOGAN  SHEAR 


• Boundary  layer  shear 

• Logan  shear  (range-  and  altitude- dependent  versions 
present  approximately  the  same  hazard) 

• Range- dependent  Kennedy  shear 

• Altitude-dependent  Kennedy  shear. 

There  are  obvious  differences  in  the  effect  any  one  shear  has  on  the  various 
aircraft  types  . We  shall,  however,  discuss  individual  aircraft  differences 
subsequently  in  Subsections  C and  D.  At  this  point  let  us  consider  one 
vehicle,  the  light  transport,  and  address  each  shear  profile  in  turn. 

A time  history  of  the  light  transport  encountering  the  boundary  layer 

shear  is  given  in  Fig.  IV- 6.  It  can  be  seen  that  the  shear  rate  (chi  /St) 

6 

is  relatively  small  even  below  200  ft  altitude  and  the  overall  effect  on 
the  aircraft  is  relatively  mild.  Note  also  that  a speed  loss  of  approxi- 
mately 8 kt  near  the  end  of  the  run  does  not  cause  any  large  angle  of  attack 
change.  This  is  due  to  the  fact  that  the  aircraft  was  in  accelerated  (less 
than  1 g)  flight  at  this  point. 

A similar  time  history  of  the  light  transport  encountering  the  range- 
dependent  Logan  shear  is  given  in  Fig.  IV- f.  Notice  that  the  aircraft 
initially  gets  high  and  1'ast  and  that  both  factors  tend  to  alleviate  the 
seriousness  of  the  shear  at  lower  altitude.  At  approximately  200  ft  (the 
altitude  at  which  the  final  abrupt  decreasing  headwind  portion  of  the  shear 
begins)  the  aircraft  is  on  the  glide  slope  but  5 kt  fast.  In  this  case 
the  aircraft  ends  up  getting  10  kt  slow  and  35  ft  below  the  glide  slope. 

It  is  also  noticeable  that  the  excursion  below  the  glide  slope  is  arrested 
chiefly  because  the  shear  rate  goes  to  zero  rather  than  because  of  the 
action  of  the  pilot  model. 

The  time  history  for  the  light  transport  encountering  the  range- dependent 
Kennedy  shear  is  shown  in  Fig.  IV-8.  The  initial  wind  encountered  by  the 
aircraft  is  a combination  of  a downdraft  and  an  increasing  headwind.  The 
effect  of  the  downdraft  is  to  drive  the  aircraft  below  the  glide  slope,  while 
the  effect  of  the  increasing  headwind  is  to  drive  the  aircraft  above  the 
glide  slope.  The  initial  result  in  this  case  is  to  drive  the  aircraft  fast 
and  below  the  glide  slope. 
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The  next  phase  of  the  shear  encounter  has  the  downdraft  turning  into 
an  updraft,  while  the  headwind  decreases  very  gradually.  The  updraft 
effect  predominates  in  this  phase  and  the  aircraft  ends  up  above  the  glide 
slope  while  being  very  near  to  its  desired  speed. 

The  final  phase  of  the  shear  encounter  finds  the  updraft  changing 
again  to  a downdraft,  and  the  headwind  rapidly  decreasing  to  a tailwind. 

Both  of  these  effects  tend  to  drive  the  aircraft  below  the  glide  slope. 

The  maximum  glide  slope  deviation  occurs  at  rouglily  200  ft  altitude,  and 
it  is  this  final  deviation  tliat  is  used  as  a measure  of  the  hazard.  The 
maximum  speed  error  is  only  16  kt  slow. 

It  should  be  noted  that  although  Ah  goes  to  a maximum  value  at  approxi- 
mately 500  ft,  tills  value  is  net  used  in  our  results.  This  is  because  at 
the  altitude  and  range  at  which  this  peak  occurs,  the  pilot  model  has  a 
lower  glide  slope  gain  per  foot  of  glide  slope  deviation  than  it  does  at 
a lower  altitude.  Hie  result  of  this  is  to  make  the  value  of  the  first 
excursion  very  dependent  upon  the  value  of  the  pilot  model  parameter  Rq 
(range  at  which  the  pilot  regulates  distance  from  glide  slope  rather  than 
angle).  Increasing  Rq  would  make  the  second  Ah  peak  the  larger  one.  Under 
these  conditions,  the  realism  of  the  first  peak  is  perhaps  questionable, 
but  the  second  peak  should  be  realistic. 

Time  histories  for  the  light  transport  encountering  the  altitude-dependent 
Logan  and  Kennedy  shears  are  shown  in  Figs.  IV-9  and  IV-10,  respectively. 

The  general  trend  of  aircraft  behavior  is  similar  to  tliat  in  the  range- 
dependent  encounters,  although  the  details  vary  due  to  the  shear /vehicle 
interaction  of  the  altitude- dependent  shears.  Thus  in  the  altitude- 
dependent  Kennedy  shear  encounter,  a 10  kt  headwind  decrease  occurs  in 
10  sec  near  500  ft;  while  in  the  range- dependent  Kennedy  shear,  the  same 
10  kt  headwind  decrease  requires  almost  20  sec.  This  higher  shear  rate  is 
typical  of  the  effects  that  occur  in  the  altitude- dependent  shears,  and  is 
the  cause  of  the  greater  deviation  from  the  glide  slope  for  all  four  air- 
craft for  these  shears. 
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C.  EFFECT  OF  VARYING  OPERATIONAL  TECHNIQUES 

After  determining  the  baseline  shear  hazard  for  the  variolas  aircraft 
types  considered,  we  then  studied  the  effects  of  varying  operational  tech- 
niques . These  operational  aspects  are  conveniently  divided  into : 

• Trim  conditions  (approach  airspeed  and  glide  slope  angle) 

• Variability  in  trim  conditions  (manual  regulation  of 
airspeed  and  flight  path) 

Both  of  these  are  explored  by  further  use  of  the  shear  hazard  metric  Ah  . 
in  conjunction  with  the  severe  Logan  and  Kennedy  shears. 

1 . Trim  Conditions 

We  begin  by  considering  the  effect  of  trimmed  approach  airspeed  because 
it  involves  the  most  fundamental  relationship  to  wind  shear  hazard  as  we 
shall  explain. 

a.  Effect  of  Trimmed  Approach  Airspeed.  To  evaluate  the  effects  of 
approach  airspeed  on  shear  hazard,  the  light  transport,  jumbo,  and  business 
jet  were  flown  through  both  range-  and  altitude-dependent  versions  of  the 
Logan  and  Kennedy  shears  at  150  kb.  The  shear  hazard  at  150  kt  was  then 
compared  to  the  shear  hazard  at  each  aircraft's  nominal  approach  airspeed. 
For  the  STOI.  commuter,  80  kt  was  chosen  as  the  higher  approach  airspeed  to 
evaluate . Figures  IV- 11  and  IV-12  show  the  results  for  the  range-  and 
altitude-dependent  versions  of  the  Logan  shear,  while  Figs.  IV- 1 3 and  IV- 1 L 
show  the  data  for  the  range-  and  altitude-dependent  versions  of  the  Kennedy 
shear . 

For  all  but  one  case  the  trend  is  that  an  airspeed  increase  decreases 
shear  hazard.  This  is  true  for  all  four  aircraft  in  the  Logan  shear,  and 
all  three  CTOL  aircraft  in  the  Kennedy  shear.  The  exception  is  the  SIDL 
commuter,  which  staffers  an  increase  in  altitude  loss  in  both  the  range- 
and  altitude-dependent  versions  of  the  Kennedy  shear.  Before  addressing 
the  anomaly,  though,  let  us  consider  the  general  trend. 

The  four  figures  illustrating  approach  airspeed  effects  show  that  Ah 

ITlcLX. 

for  all  the  aircraft  generally  falls  along  an  inverse  speed  trend  line.  To 
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explain  this  phenomenon  we  can  utilize  a feature  of  linear  analysis . Recall 

from  the  previous  section  that  the  stability  derivatives  which  relate 

incremental  normal  acceleration  to  velocity  increments  of  tailwind  (u  ) 

IS 

and  downdraft  (w  ) are  Z and  Z , respectively.  Considering  the  sensitivity 
g u w 

to  horizontal  gust  components,  we  should  expect  to  see  that  the  size  of  an 
altitude  change  corresponds  to  the  magnitude  of  Z^.  This  is  precisely  the 
trend  we  obtain  since 


i.e.,  Z is  inversely  proportional  to  airspeed  as  is  the  resulting  altitude 

excursion.  Ah  for  nearly  all  cases  simulated. 

’ max 

The  correlation  of  altitude  excursion  with  airspeed  is  an  important 
result  for  it  establishes  the  basic  first  order  effect  of  a given  wind 
shear  on  nearly  any  conventional  airplane  (which  includes  the  STOL  commuter). 

The  counter  example  to  the  general  effect  of  airspeed  on  shear  hazard 
is  the  STOL  commuter  in  the  Kennedy  shear  (Figs.  IV-13  and  IV-lL).  The 
source  of  the  anomaly  can  be  traced  first  to  the  vertical  gust  upset  which 
was  present  in  the  Kennedy  shear  and  not  the  Logan  shear.  In  addition, 
there  is  an  increased  sensitivity  of  the  STOL  commuter  to  vertical  gusts 
which  is  related  to  a pilot- vehicle  interaction. 

The  STOL  commuter  pilot-vehicle  combination  is  inherently  more  sensitive 
to  vertical  gusts  than  the  CTOL  examples  because  of  a strong  flight  path- 
to-airspeed  cross  coupling  quality  which  produces  a low  damping  ratio  oscil- 
latory mode  when  the  flight  path  loop  is  closed.  This  difference  is  directly 
related  to  the  effective  thrust  angle  which  is  U5  deg  for  the  STOL  commuter 
and  nearly  zero  for  the  CTOL  examples.  Figure  IV-15  shows  how  the  cross 

coupling  effect  manifests  itself  when  excited  by  a step  w . The  altitude 

6 

excursion  is  normalized  with  the  closed  loop  flight  path  natural  frequency 
(which  reflects  the  precision  of  the  flight  path  loop  closure)  and  heave 
damping,  Z , . Thus,  the  normalized  altitude  excursion  in  Fig.  IV-15  reflects 
a direct  comparison  of  flight  path-to-airspeed  cross  coupling  and  explains 
the  high  sensitivity  of  the  STOL  commuter  to  vertical  gusts. 
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The  reason  for  the  STOL's  increased  shear  hazard  at  80  kt  is  attribu- 
table to  the  increased  heave  damping,  Z^,  which  is  proportional  to  airspeed. 

To  sum  up  the  results,  then,  the  major  contributor  to  wind  shear  hazard 
is  the  horizontal  gust  sensitivity,  Z , which  is  inversely  proportional  to 
airspeed.  This  explains  the  gross  difference  in  hazard  between  the  slower 
STOL  commuter  and  the  faster  CTOL  aircraft  and  the  difference  for  increased 
approach  speeds  in  the  CTOL  aircraft  alone.  However,  where  a significant 
vertical  gust  component  is  involved,  the  STOL  commuter  does  not  necessarily 
enjoy  the  reduction  in  hazard  due  to  increased  airspeed. 

b . Effects  of  Glide  Slope  Angle  on  the  STOL  Commuter  Aircraft . The 

effects  of  glide  slope  angle  on  shear  hazard  were  tested  for  the  STOL 

commuter  aircraft  using  a three  and  six  deg  glide  slope.  The  aircraft  model 

was  slightly  modified  for  this  series  of  tests . As  we  shall  discuss  shortly, 

the  thrust  offset  of  the  STOL  commuter  appeared  unrealistically  large  and 

resulted  in  a large  adverse  effect  on  Ah  .In  order  to  prevent  thrust 
° max 

offset  from  dominating  the  study  of  glide  slope  angle,  the  offset  was 
nulled  (by  setting  the  thrust  line  coincident  with  the  eg). 

The  data  for  the  Logan  shear,  both  range-  and  altitude-dependent,  are 
shown  in  Fig.  IV- lG  and  the  data  for  the  Kennedy  shear,  both  range-  and 
altitude- dependent,  in  Fig.  IV-17.  The  results  do  not  show  a consistent 
trend.  In  the  Logan  shear,  use  of  a 5 deg  glide  slope  reduces  the  shear 
hazard  by  a factor  of  one-half.  In  the  Kennedy  shear,  the  effect  of  glide 
slope  variation  is  small  compared  to  the  effect  of  shear  dependency. 

To  explain  these  results  it  is  useful  to  refer  to  the  closed  loop 
frequency  response  plots  from  Section  III  (Fig.  Ill-  6).  The  STOL  commuter 
Ah/u  frequency  response  is  fairly  sharply  tuned  to  a frequency  of  approxi- 
mately  0.5  rad/sec.  For  a shear  cycle  which  occurs  over  a height  of,  say, 

500  ft  and  an  average  sink  rate  of  10  ft/sec,  the  dominant  spectral  power 
would  be  centered  at  about  2 n * or  0.12  rad/sec.  Therefore,  if  we 
were  to  approach  at  one-half  the  sink  rate,  or  5 ft/sec,  the  dominant 
spectral  power  of  the  disturbance  would  shift  downward  to  0.06  rad/sec. 

The  slope  of  the  Ah/u^  frequency  response  at  low  frequencies  would  predict 
an  attenuation  of  Ah  equal  to  the  spectral  shift.  This  is  precisely  the 
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Commuter  in  Logan  Shear 


case  in  the  Logan  shear  in  which  there  is  no  vertical  gust  component.  The 
variation  from  6 deg  glide  slope  angle  to  3 deg  reduces  the  abruptness  of 
the  wind  shear  disturbance  and  thus  the  pilot-vehicle  system  more  effectively 
filters  out  the  shear  effect,  in  fact,  by  a factor  of  one- half. 

In  the  Kennedy  shear,  the  results  are  not  so  easily  explained  because 

of  the  presence  of  a significant  vertical  gust  component.  While  there  is  a 

spectral  shift  in  the  disturbance  which  corresponds  to  varying  sink  rate, 

there  is  an  interaction  of  vertical  and  horizontal  gust  components  which  is 

offsetting.  Although  the  horizontal  disturbance  component  is  more  effectively 

filtered  out  on  the  3 deg  glide  slope,  the  vertical  component  remains  just 

as  effective  as  for  the  6 deg  according  to  the  Ah/w  frequency  response 

g 

plot  in  Fig.  Ill- 7 which  is  flat  at  low  frequencies. 

2.  Variability  in  Trim  Conditions 

a.  Effects  of  Airspeed  Regulation.  The  characteristics  of  the  airspeed 
regulation  loop  were  evaluated  during  the  simulation  by  altering  both  the 
airspeed  gain,  K , and  the  airspeed  indifference  threshold,  (=UH)*  I* 
was  suspected  that  allowing  the  pilot  model  to  track  airspeed  more  tightly 
or  loosely  could  have  an  impact  on  shear  hazard.  The  nominal  5 kt  airspeed 
indifference  threshold  was  varied  from  0 to  10  kt,  and  the  airspeed  loop 
gain,  Ku,  was  -varied  from  one-half  to  twice  its  nominal  value.  However, 
the  airspeed  loop  gain  variation  is  represented  in  terms  of  airspeed  loop 
crossover  frequency,  toCu,  for  more  generality. 

The  results  of  varying  airspeed  loop  threshold  effects  on  shear  hazard 
are  presented  in  Figs.  IV-18  and  IV-19  for  the  light  transport,  jumbo,  and 
STOL  commuter.  The  light  transport  and  jumbo  aircraft  show  a decrease  in 
shear  hazard  as  the  airspeed  indifference  threshold  is  decreased  for  both 
versions  of  the  Kennedy  shear.  The  STOL  commuter  aircraft  shows  a similar 
trend  for  the  range-dependent  Kennedy  shear,  but  for  the  altitude-dependent 
Kennedy  shear  there  is  no  significant  change  in  shear  hazard  for  a decrease 
in  airspeed  indifference  threshold. 

The  data  reflecting  the  effects  of  airspeed  loop  crossover  frequency 
(loop  gain)  are  shown  in  Figs.  IV-20  and  IV-21  for  the  range-  and  altitude- 
dependent  versions  of  the  Kennedy  shear,  pata  are  shown  only  for  the  light 
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-igure  IV-20.  Effects  of  Airspeed  Loop  Crossover  Frequency  in  Kennedy  Shear 

( Range- Dependent ) 


transport  and  jumbo  aircraft.  The  STOL  commuter  airspeed  loop  dynamics 
would  not  permit  a simple  doubling  or  halving  of  the  airspeed  loop  gain 
without  a departure  from  the  loop  closure  rules  adopted,  therefore  it  was 
dropped  in  this  case.  The  plotted  data  in  Figs.  IV-20  and  IV-21  show  that 
for  the  light  transport  and  jumbo  aircraft  the  precision  of  speed  regulation 
has  a strong  effect  on  shear  hazard.  With  tighter  control  of  airspeed 
during  the  shear  encounter,  less  lift  is  lost,  and  the  subsequent  glide 
slope  deviation  is  reduced. 

The  data  for  the  effects  of  airspeed  threshold  and  airspeed  loop  gain 
are  consistent  in  that  they  reflect  the  fact  that  an  increase  in  the  precision 
of  airspeed  regulation  results  in  a reduction  of  shear  hazard  for  the  light 
transport  and  jumbo  aircraft.  It  appears  that  either  decreasing  the 
airspeed  indifference  threshold  or  increasing  the  airspeed  loop  gain  are 
roughly  equivalent  ways  of  increasing  airspeed  regulation  although  both 
imply  markedly  increased  concentration  of  attention  (and  devotion  of  work- 
load) to  airspeed  regulation  by  a human  pilot. 

b.  Effects  of  Flight  Path  Regulation.  Since  the  shear  hazard  metric 
used  in  this  study  was  defined  to  be  aircraft  deviation  from  the  glide 
slope,  it  is  reasonable  to  expect  the  shear  hazard  to  vary  as  a function  of 
how  tightly  the  pilot  closes  a flight  path  loop.  During  the  simulation, 
flight  path  loop  tightness  was  adjusted  by  changing  the  value  of  the 
parameter  in  the  flight  path  loop.  The  light  transport,  jumbo,  and  STOL 
commuter  were  flown  through  range-  and  altitude- dependent  versions  of  the 
Kennedy  shear  with  the  values  of  alternatively  set  to  half  and  twice  the 
nominal  values . While  K^  was  the  variable  adjusted,  it  is  more  generally 
useful  to  present  the  data  in  terms  of  flight  path  loop  crossover  frequency, 
a>C(i  (i.e.,  closed  loop  system  bandwidth).  The  data  for  the  range-  and 
altitude- dependent  Kennedy  shear  is  shown  in  Figs.  IV-22  and  IV-23,  res- 
pectively. The  higher  crossover  frequency  of  0.5  rad/sec  represents  a 
workload  increase  by  a human  pilot  that,  while  high,  is  certainly  attain- 
able, especially  in  later  stages  of  the  approach.  A crossover  frequency  of 
0.15  rad/sec  might  be  used  in  early  stages  of  the  approach  prior  to  a 
sensitive  glide  slope  indication. 
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The  shear  model  dependency  appears  to  have  a strong  influence  on  the 
overall  effectiveness  of  flight  path  regulation  in  reducing  shear  hazard. 

For  the  range- dependent  Kennedy  shear  the  Ah^  excursions  for  each  air- 
craft are  similar  in  magnitude  and  approximately  inversely  proportional  to 
flight  path  crossover  frequency,  u)  . The  altitude-dependent  shear,  however, 
shows  a wider  disparity  of  hazard  among  aircraft  types  and  less  variation 
with  flight  path  loop  tightness.  This  is  an  example  of  how  greatly  the 
feedback  effect  of  an  altitude- dependent  shear  (as  shown  in  Fig.  IV- 3) 
can  modify  simulator  results.  Notice  that  the  altitude-dependent  shear 
does  not  always  result  in  a worst  case  as  shown  by  the  light  transport  with 
a loose  flight  path  loop  (i.e.,  a low  crossover  frequency,  a>c^) . 

D.  EFFECT  OF  VARYING  PROPULSION  SYSTEM  FEATURES 

A few  propulsion  system  features  were  studied  using  the  wind  shear 
simulation.  These  features  were  mainly  performance-related  and  therefore 
strongly  tied  to  propulsion  system  dynamics.  The  parameters  investigated 
include: 


• Thrust  available 

• Thrust  lag 

• Thrust  offset. 

1 . Thrust  Available 

Maximum  thrust  available  is  normally  considered  to  be  the  basic  ingre- 
dient in  climb  performance,  thus,  it  was  natural  to  explore  its  effect  with 
respect  to  a wind  shear  encounter. 

The  means  of  exploring  the  effect  of  thrust  available  was  simply  to 
limit  it  during  a simulated  shear  encounter.  This  was  done  using  the  light 
transport  in  the  altitude-dependent  Kennedy  shear.  The  trim  thrust- to-weight 
ratio  was  0.1 5 and  the  maximum  incremental  thrust- to-weight  demanded  by  the 
mathematical  p. lot  model  was  + 0.1 4.  The  wind  shear  hazard,  Ahmay,  is 
plotted  in  Fig.  IV-2U  for  available  thrust  levels  below  the  maximum  demanded. 
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The  results  of  this  investigation  show  that  limited  thrust-to-weight 

ratio  has  a minor  effect  in  contributing  to  an  increase  in  wind  shear 

hazard  in  the  case  studied.  The  reason  for  this  is  that  the  CTOL  aircraft 

depend  little  upon  thrust  for  short-term  flight  path  response.  Instead, 

thrust  is  modulated  in  a longer  term  context  to  preserve  adequate  trimmed 

airspeed  and  thus  to  insure  effective  short  term  flight  path  response  via 

attitude  control  acting  through  the  airspeed- heave  damping  product  U Z 

o w 

discussed  previously  in  Section  II. 

2.  Thrust  Lag 

The  lag  between  a throttle  input  and  thrust  output  was  varied  T >m  the 
nominal  value  of  each  aircraft  to  a common  value  of  0.5  sec. 

The  results  are  shown  in  Figs.  IV-25  and  IV-26  for  the  range-  and 
altitude-dependent  Logan  shears  and  Figs.  IV -27  and  IV-28  for  the  range- 
and  altitude-dependent  Kennedy  shears,  respectively.  The  thrust  lag  reduc- 
tion (or  increase  in  the  case  of  the  light  transport)  from  its  nominal 
value  has  a minor  effect  in  contributing  to  a reduction  (or  increase)  in 
wind  shear  hazard  for  the  CTOL  aircraft.  Again,  as  in  the  case  of  available 
thrust-to-weight  ratio,  the  reason  is  that  CTOL  aircraft  depend  little 
upon  thrust  for  short-term  flight  path  response.  The  jumbo  shows  the  most 
exaggerated  effect  on  hazard,  because  nominal  lag  is  5 sec. 

In  the  case  of  the  STOL  commuter,  the  thrust  lag  is  more  directly 
involved  in  AhmRy,  since  the  flight  path  is  primarily  controlled  by  the 
throttle  in  the  STOL  technique  examined  here. 

5.  Thrust  Offset 

There  are  two  pi telling  moment  effects  related  to  vertical  displacement 
of  the  effective  thrust  line.  The  first,  pitching  response  to  thrust,  can 
be  quantified  as  the  ratio  of  pitching  acceleration  to  longitudinal  accelera- 
tion for  a throttle  input,  M^/X^.  The  second  effect,  pitching  response 
to  airspeed,  can  be  quantified  by  the  stability  derivative  M^. 

The  sources  of  Mf)T/X&T  and  Mu  are  illustrated  in  Fig.  IV-29.  A jet 
engine  is  depicted,  however,  a propeller-driven  propulsion  system  shares 
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Light  Transport 


igure  17-26.  Effects  of  Thrust  Lag  in  Logan  Shear 
( Altitude-Dependent ) 


Light  Transport 
Business  Jet 


Figure  IV-28.  Effects  of  Thrust  Lag  in  Kennedy  Shear 
(Altitude -Dependent ) 
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Figure  IV-29.  Contributions  of  Thrust  Offset  to  Pitching  Moment  Derivatives 
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the  same  features.  Note  that  a downward  displacement  of  the  thrust  line 
from  the  vehicle  eg  produces  the  following  results: 

• An  increase  in  thrust  produces  a pitch-up  moment 
from  direct  thrust  effect 

• An  increase  in  tailwind  produces  a pitch-up  moment 
from  decreased  ram  drag. 

Both  of  these  effects  will  be  countered  somewhat  by  the  pilot's  pitch 
attitude  loop  which  tries  to  hold  attitude  constant  in  the  absence  of  pilot 
commanded  pitch  inputs . However,  there  is  sufficient  vehicle  response  in 
the  short  term  for  the  and  M&T/X&T  effects  to  be  noticeable.  As  the 
tightness  of  the  attitude  loop  increases,  these  effects  will  become  less 
noticeable . 

Variations  in  thrust  offset  were  evaluated  during  the  simulation  because 
there  was  some  doubt  about  the  validity  of  the  thrust  offsets  for  the  light 
transport  and  STOL  commuter.  Comments  from  pilots  who  have  flown  both  the 
actual  aircraft  and  simulations  with  these  nominal  engine  offsets  indicate 
that  the  effective  thrust  offset  for  the  simulated  vehicles  was  excessive. 
The  effective  offset  is  difficult  to  estimate  and  can  be  considerably 
different  from  the  geometric  offset  due  to  engine  flow  and  downwash  effects 
on  the  horizontal  tail.  An  additional  reason  for  varying  thrust  offset  is 
that  a pilot  could  largely  cancel  the  thrust  offset  effects  by  using  a 
tighter  pitch  attitude  loop  or,  more  likely,  by  establishing  a throttle-to- 
elevator  crossfeed.  Table  IV-1  lists  the  thrust  offsets  and  the  Mg^/Xg^, 
ratio  for  each  aircraft. 

The  effects  of  thrust  offset  on  shear  hazard  were  evaluated  by  varying 
the  thrust  offset  from  its  nominal  value  to  zero.  The  results  for  the 
range-  and  altitude-dependent  versions  of  the  Logan  and  Kennedy  shears 
are  given  in  Figs.  IV-30  through  IV- 33,  respectively. 

The  trend  of  the  data  is  consistent  in  that  increasing  the  pitch-to- 
thrust  coupling  by  displacing  the  thrust  line  downward  decreases  the  shear 
hazard  by  producing  an  immediate  favorable  pitch-up  response  to  thrust 
increase  and  tailwind  increase. 
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TABLE  IV- 1 

THHIEI  OFFSET  AMD  PITCffiK  1®^  ro  THEUSI 


AIRCRAFT 

Jumbo 

Light  Transport 
Business  Jet 
STOL  Commuter 


OFFSET 

(Positive  Below  cr)  ft 


10 

5 

-1.5 

-2.9 


M6t/x5t  rad/ft 

.OO55 
.0169 
-.0318 
-.04 27 
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Jumbo 


E.  SUMMARY 


A number  of  features  have  been  examined  regarding  wind  shear  effects 
in  the  approach  flight  phase.  Let  us  summarize  the  results  in  order  to 
separate  the  first  order  effects  from  the  higher  order  effects. 

First,  our  linear  analysis  helped  to  point  out  the  important  pilot- 

vehicle  features.  Recall  that  the  altitude  excursion  transfer  relationships 

showed  that  the  horizontal  gust  components  are  most  effective  in  the  spectral 

region  bounded  from  below  by  airspeed  regulation  and  from  above  by  flight 

path  regulation  bandwidth  (l/T  and  u>d).  The  magnitude  of  the  closed  loop 

response  is  inversely  proportional  to  airspeed  itself  and  the  frequency  of 

the  flight  path  regulation  mode  cu^.  Vertical  gust  components  produce  a 

response  from  very  low  frequencies  to  the  frequency  co..  The  amplitude  of 

a p 

the  response  is  proportional  to  heave  damping,  Z , and  the  inverse  of  ok  . 

w d 

The  Logan  wind  shear  helps  to  illustrate  the  effect  of  a large  amplitude, 

strictly  horizontal  gust  profile  on  a wide  variety  of  passenger  aircraft. 

Figures  IV-34  and  IV- 35  show  the  hazard  metric  Ah  for  both  the  range- 

max  ° 

and  altitude-dependent  versions  of  the  shear.  Key  propulsion  features  have 

been  made  common  among  the  aircraft  to  minimize  second  order  effects,  i.e., 
thrust  offset  is  zero  and  thrust  lag  is  0.5  sec  for  al 1 aircraft.  The  data 
follow  the  inverse  airspeed  trend  line  well.  In  addition,  the  hazard  metrics 
for  the  CTOL  group  of  aircraft  coalesce  even  better  than  the  results  in 
previous  Figs.  IV-11  through  IV- 1 4 demonstrated.  This  illustrates  the 
limited  importance  of  increased  airspeed,  whether  viewed  as  a design  feature 
or  an  operational  technique,  in  reducing  the  shear  hazard. 
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Light  Transport 


Summary  of  Effects  of  Airspeed  in  the  Logan  Shear 
( Range-Dependent ) 
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(Altitude-Dependent ) 


SECTION  V 


TAKEOFF  RESULTS 


A large  number  of  takeoff  cases  were  tested.  The  range-  and  altitude- 
dependent  shears  were  tested  for  each  of  the  four  aircraft  with  variations 
in  trim  speed,  initial  rate  of  climb,  and  parameters  in  the  pilot  model. 
This  section  is  divide  d into  three  subsections  each  of  which  is  concerned 
with  a specific  topic.  The  first  subsection  will  discuss  the  takeoff  shear 
hazard  criteria.  Subsection  B will  discuss  the  shear  effects  for  the  base- 
line vehicles.  Subsection  C will  discuss  potential  operational  techniques 
which  have  an  impact  on  shear  hazard. 


A.  SHEAR  HAZARD  CRITERIA 

The  shear  hazard  metric  used  for  the  approach  cases,  maximum  deviation 
from  the  glide  slope,  is  clearly  not  suitable  for  a takeoff  shear  hazard 
criteria.  During  takeoff  there  is  not  a glide  slope  for  the  airplane  to 
track.  Rather,  the  primary  tracking  task  involves  airspeed  control  using 
pitch  attitude.  There  is  no  throttle  control  loop  since  the  aircraft  is 
assumed  to  remain  at  takeoff  power.  The  takeoff  shear  hazard  metric  decided 
upon  was  the  minimum  rate  of  climb  which  resulted  from  the  shear  encounter. 
Another  possible  metric  considered  was  the  minimum  speed  resulting  from 
the  shear  encounter.  However,  the  aircraft  is  not  in  a trimmed  state  (i.e., 
1 g flight)  during  the  encounter,  and  the  trim  angle  of  at+" 'k  does  not 
approach  a limiting  case  (stall).  Thus  minimum  speed  during  the  shear 
encounter  was  not  considered  to  be  a critical  hazard  metric. 


As  indicated  in  Section  II,  the  procedure  used  during  the  takeoff  tests 
was  to  trim  the  aircraft  at  a height  of  30  ft.  The  aircraft  airspeed  was  set 
at  either  a published  Vg  speed  or  at  1.3  vs^a]_]_  for  the  takeoff  configuration. 
This  shear  was  oriented  so  that  the  aircraft  was  taking  off  into  a headwind. 
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The  takeoff  model  utilizes  an  airspeed-to-pitch  attitude  control  feed- 
back with  a 3 kt  airspeed  threshold.  The  power  remains  constant  at  the 
takeoff  power  setting.  The  rules  for  adjusting  the  pitch  attitude  and 
airspeed  loops  are,  as  indicated  in  Section  II,  identical  to  those  for  the 
approach  case. 

A brief  restatement  of  the  problem  encountered  by  an  aircraft  taking 
off  in  the  presence  of  a decreasing  headwind  shear  will  help  to  illuminate 
the  results.  Typically,  the  takeoff  aircraft  is  at  or  near  its  best  rate 
of  climb  speed  (V  ) at  a 30  ft  altitude  and,  following  a loss  in  speed, 

y 

the  airplane  must  reduce  the  pitch  attitude  (and  angle  of  attack)  to  accel- 
erate back  to  V . Any  attempt  to  fly  at  a speed  less  than  V would  result 

y y 

in  a decreased  rate  of  climb.  Additionally,  any  attempt  to  increase  the 
rate  of  climb  by  increasing  pitch  attitude  would  only  cause  the  aircraft 
to  decelerate  and  would  ultimately  reduce  the  rate  of  climb.  It  should 
also  be  noted  that  the  immediate  problem  of  an  aircraft  upon  encountering 
the  decreasing  headwind  shear  is  not  proximity  to  stall  but  rather  the  loss 
of  rate  of  climb  due  to  the  airspeed  loss. 

B.  SHEAR  EFFECTS  FOR  THE  BASELINE  VEHICLES 

As  indicated  in  the  previous  section,  the  shear  hazard  metric  used  during 
the  takeoff  simulations  was  the  minimum  rate  of  climb  encountered  during 
the  shear  encounter.  Figures  V-l  and  V-2  depict  the  results  for  the  range- 
and  altitude-dependent  Logan  and  Kennedy  shears  for  the  baseline  aircraft. 
Only  the  severest  shear,  the  Kennedy  shear,  was  tested  among  all  aircraft, 
because  of  time  limitations  during  the  simulation.  It  can  be  seen  from 
the  data  that,  just  as  for  the  approach  cases,  we  can  rank  the  shears  in 
terms  of  increasing  severity: 

( 1 ) Boundary  layer 

( 2 ) Logan 

( 3 ) Kennedy 

It  can  also  be  seen  that,  of  the  three  shears,  the  only  shear  capable  of 
reducing  the  aircraft  rate  of  climb  to  zero  is  the  Kennedy  shear.  Table  V-l 
lists  the  nominal  trimmed  takeoff  rate  of  climb  for  each  of  the  aircraft. 
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Takeoff  Results  in  Altitude-Dependent  Shears 


TABLE  V-1 


NOMINAL  AIRCRAFT  TRIM  STATES 


AIRCRAFT 

TAKEOFF  AIRSPEED 

TRIM  h 

(kt) 

(ft/sec 

Jumbo 

187 

38.7 

Light  Transport 

l6l 

37-4 

Business  Jet 

135 

62.5 

STOL  Commuter 

78 

21.5 

It  can  be  seen  that  the  relatively  good  performance  of  the  business  jet  in 
the  range-dependent  Kennedy  shear  is  due  to  its  large  initial  trimmed  rate 
of  climb.  In  fact,  each  of  the  aircraft  loses  approximately  30  ft/sec 
rate  of  climb  during  an  encounter  with  the  range-dependent  Kennedy  shear. 

It  should  be  noted  that  the  alti tilde- dependent  shears  have  a shear-pilot- 
vehicle  interaction  on  takeoff  which  is  similar  to  the  approach  interaction 
discussed  earlier.  The  effect  of  this  interaction  is  to  give  most  of  the 
aircraft  virtually  the  same  shear  hazard  for  the  altitude- dependent  case 
regardless  of  their  initial  rate  of  climb.  Comparing  the  business  jet 
shear  hazard  in  the  range-  and  altitude- dependent  Kennedy  shear  cases,  we 
see  that  this  aircraft  suffers  substantially  more  for  the  altitude- dependent 
case.  This  is  due  to  the  fact  that  this  aircraft's  higher  initial  rate-of- 
climb  in  an  altitude-dependent  shear  results  in  a more  rapid  rate  of  shear 
encounter  than  for  the  other  aircraft.  Thus  it  experiences  a more  severe 
shear. 

The  main  conclusions  for  the  shear  effect  on  the  baseline  aircraft 
are  that  the  shears  maintain  their  same  order  of  severity  in  the  takeoff 
cases  as  they  do  in  the  approach  cases,  and  that  a shear- vehicle- pilot 
interaction  exists  for  the  altitude- dependent  Kennedy  shear  which  tends  to 
create  a uniform  shear  hazard  for  the  aircraft  tested. 

C.  POTENTIAL  OPERATIONAL  TECHNIQUES 

The  number  of  potential  operating  techniques  which  could  be  employed  to 
reduce  takeoff  shear  hazard  is  limited.  This  is  because  the  aircraft  is 
operating  at  or  near  its  airspeed  for  best  rate  of  climb  during  takeoff. 
Additionally,  the  takeoff  pilot  model  has  only  an  airspeed- to- pitch  attitude 
control  feedback  loop,  and  the  possible  alternative  piloting  techniques 
available  are  limite’ . During  the  simulation,  al  1 pilot  parameters  as 
well  as  the  aircraft  operating  conditions  were  varied  to  determine  their 
effect  on  shear  hazard.  The  parameters  determined  to  be  of  importance  were: 

( 1 ) Takeoff  airspeed 

(2)  Airspeed  regulation  characteristics. 

These  will  be  discussed  in  order. 


I 
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Clearly,  decreasing  the  takeoff  speed  is  not  a real  is tic  option  for  a 
pilot  prior  to  a shear  encounter . It  was  found  that  increasing  the  takeoff 
speed  also  has  little,  if  any,  beneficial  effects  during  the  encounter.  This 
is  shown  in  Figs . V-3  and  V-4  which  depict  the  Kennedy  shear  hazard  for  the 
business  jet,  light  transport,  and  STOb  commuter  as  a function  of  airspeed. 
This  data  indicates  that  increasing  the  takeoff  airspeed  will  cause  a range- 
dependent  shear  to  be  encountered  more  quickly.  likewise,  for  an  altitude- 
dependent  shear,  the  subsequent  decrease  in  initial  rate  of  climb  which  is 
a consequence  of  the  increased  takeoff  airspeed  will  result  in  having  less 
rate  of  climb  margin  available  during  the  shear  encounter. 

The  only  other  potential  operational  technique  we  examined  was  airspeed 
regulation.  This  showed  a significant  effect  on  shear  hazard  but  the  effect 
is  opposite  to  that  seen  in  the  landing,  because  the  hazard  metric  is  dif- 
ferent. Figures  V-5  and  V-6  show  the  shear  hazard  for  the  light  transport, 
jumbo,  and  STOL  commuter  as  a function  of  airspeed  loop  crossover  frequency 
(closed  loop  bandwidth).  As  in  the  approach  section,  the  airspeed  loop 
gain  was  varied  by  a factor  of  2 to  determine  its  effect,  however,  the  data 
is  presented  in  terms  of  the  airspeed  loop  crossover  frequency  for  a more 
general  interpretation.  We  can  see  that  increasing  the  tightness  of  the 
airspeed  loop  (i.e.,  increasing  the  crossover  frequency)  increases  the 
takeoff  shear  hazard  whereas  decreasing  the  airspeed  loop  crossover  fre- 
quency decreases  the  takeoff  shear  hazard. 

This  effect  is  a result  of  the  actions  of  the  pilot  when  encountering  a 
decreasing  headwind  shear.  As  the  airspeed  decreases,  the  pilot  will  lower 
the  nose  of  the  aircraft  to  regain  airspeed.  Thus,  not  only  does  the 
action  of  the  shear  tend  to  reduce  the  aircraft  rate  of  climb,  but  also 
the  action  of  the  pilot  model  in  tracking  airspeed  tends  to  also  reduce  the 
rate  of  climb.  By  increasing  the  airspeed  loop  crossover  frequency  (in- 
creasing the  loop  tightness),  a larger  nose-down  pitch  change  for  each 
increment  of  airspeed  loss  will  result.  This  larger  pitch-down  change  will 
result  in  a larger  loss  in  rate  of  climb.  Likewise,  by  decreasing  the 
airspeed  loop  crossover  frequency  (decreasing  the  loop  tightness)  there 
will  be  less  of  a pitch  change  during  the  shear  encounter  and  thus  the 
pilot  model  will  have  less  of  an  impact  on  reducing  the  aircraft's  rate  of 
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Effects  of  Takeoff  Airspeed  in  Kennedy  Shear 


Effects  of  Takeoff  Airspeed  in  Kennedy 
( Altitude- Dependent ' 


Figure  V-6.  Effects  of  Airspeed  Loop  Crossover  Frequency  in  Kennedy  Shear 

(Altitude-Dependent ) 


climb.  The  implication  here,  clearly,  is  that  a definite  tradeoff  exists 
between  tracking  airspeed  loosely  and  suffering  a loss  in  rate  of 
because  airspeed  is  tracked  tightly.  It  appears  that  to  minimize  the  loss 
in  rate  of  climb  during  takeoff  the  airspeed  variability  must  be  allowed 
to  be  larger  than  would  otherwise  be  tolerated. 


SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  section  will  be  divided  into  three  subsections:  approach  conclu- 
sions, takeoff  conclusions,  and  recommendations. 

A.  APPROACH  CONCLUSIONS 

The  following  summary  of  conclusions  for  the  approach  flight  phase  is 
divided  into  four  major  areas,  (i)  validity  and  limitations  of  the  simu- 
lation, (ii)  results  of  linear  analysis,  (iii)  results  from  studies  of  base- 
line aircraft  conf igurations , * and  (iv)  pilot-vehicle  features  which  affect 
the  shear  hazard. 

1 . Simulation  Validity  and  limitations 

• An  unmanned  simulation  utilizing  a suitable  pilot  model  can  be  an 
effective  tool  for  investigating  pilot- vehicle- shear  interactions. 

• One  simple  and  easily  quantifiable  measure  of  the  hazard  to  an 
aircraft  on  approach  is  its  maximum  vertical  deviation  from  the 
glide  slope  due  to  wind  disturbances.  This  hazard  metric  was  used 
throughout  the  approach  study. 

• Other  parameters  such  as  speed  and  angle  of  attack  were  considered 
as  shear  hazard  metrics  but  were  discarded,  either  because  of  in- 
sensitivity to  variations  of  aircraft  parameters  or  because  they 
did  not  approach  limiting  values  as  did  altitude  deviation. 

• Shear  model  characteristics  can  couple  with  vehicle  dynamics  to 
produce  what  may  be  exaggerated  results.  The  effects  of  shear 
dependency  (altitude,  range,  time)  must  be  recognized  and 
appreciated. 

* The  baseline  aircraft  configurations  are  the  unaltered  models  of  the  four 
aircraft  involved  in  this  study. 
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2.  Results  of  Linear  Analysis 

• Linear  analysis  was  a powerful  tool  far  (i)  reducing  the  complex 
simulation  models  used  in  this  study  to  a minimum  number  of 
essential  parameters,  (ii)  showing  the  specific  interrelation  of 
those  parameters  to  wind  shear  hazard,  and  (iii)  showing  some 
effects  of  shear  dependency  on  altitude,  range,  and  time. 

• For  horizontal  gust  (or  wind  shear)  components  linear  analysis 
predicts  that  the  pilot- vehicle- shear  sensitivity  is  proportional 
to: 

(i)  Inverse  airspeed 

(ii)  Inverse  airspeed  loop  crossover  frequency 
(loop  tightness) 

(iii)  Inverse  flight  path  loop  crossover  frequency  squared 

• For  vertical  gust  components  linear  analysis  predicts  that  the 
pilot-vehicle-shear  sensitivity  is  proportional  to; 

(i)  Airspeed 

(ii)  Lift  curve  slope  (c^) 

(iii)  Inverse  wing  loading 

(iv)  Inverse  flight  path  loop  crossover  frequency  squared 

• The  spectral  regions  of  interest  for  wind  shear  hazard  are 
below  the  flight  path  loop  crossover  frequency. 

3.  Results  From  Studies  of  Baseline  Aircraft  Configurations 

• The  boundary  layer  shear  posed  the  least  hazard  to  all  aircraft. 
This  is  because  it  had  very  low  shear  rates  at  altitudes  above 
200  ft,  and  below  200  ft  the  higher  shear  rates  did  not  have 
sufficient  time  to  create  a large  altitude  deviation. 

• The  more  severe  Logan  shear  presents  a substantially  greater 
hazard  to  the  STOL  commuter  than  to  the  other  three  aircraft. 

This  is  because  of  the  STOL  commuter's  lower  approach  speed  and 
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its  consequently  higher  normal  acceleration  increment  due  to 
airspeed  change  (i.e.,  the  stability  derivative  Zu  is  greater). 

• For  the  most  severe  Kennedy  shear,  the  altitude  deviations  vary 
over  a + 15  percent  band  with  an  average  of  65  ft.  In  this  case 
the  effects  of  horizontal  shears  and  vertical  gusts  combine  in 
a complex  way  such  that  no  single  key  factor  can  be  identified 
for  determining  shear  hazard. 


4.  Reduction  in  Shear  Hazard  Due  to  Pilot-Vehicle  Features 


• The  methods  used  to  reduce  shear  hazard  can  be  lumped  into  three 
categories,  (i)  increasing  approach  speed,  (ii)  increasing  pilot 
tracking  precision  on  flight  path  and/or  airspeed,  and  (iii) 
altering  the  propulsion  system  dynamics. 


• Increasing  approach  speed  has  two  main  effects  for  range- 
dependent  horizontal  shears 


(i)  Decreasing  the  normal  acceleration  increment  due 
to  airspeed  change  (Z  ) 

(ii)  Increasing  the  shear  rate  (du g/^t) 


The  first  effect  causes  a decrease  in  the  vehicle  shear  response 
while  the  second  causes  an  increase,  i.e.,  the  effects  cancel  to 
some  extent.  The  overall  effect  of  increasing  airspeed  for  the 
Logan  and  Kennedy  shears,  however,  was  to  generally  reduce  the 
shear  hazard. 


• Increasing  approach  speed  has  two  main  effects  for  range- 
dependent  vertical  shears 

(i)  Increasing  the  normal  acceleration  increment 

due  to  vertical  speed  change  (i.e.,  heave  damping 

°rZ„) 

(ii)  Decreasing  the  time  the  aircraft  spends  in  the 
vertical  shear 


For  the  Kennedy  shear  these  effects  appeared  to  cancel. 
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• For  roughly  comparable  increases  in  precision,  better  tracking 
of  flight  path  results  in  a greater  shear  hazard  reduction 

(30  percent  reduction  in  altitude  deviation  for  the  light 
transport  in  the  Kennedy  shear)  than  does  more  precise  tracking 
of  airspeed  (18  percent  reduction  for  the  same  case).  This  is 
because  (i)  airspeed  has  an  indirect  effect  on  altitude  deviation 
while  flight  path  has  a direct  effect  on  it,  and  (ii)  the 
presence  of  propulsion  system  lag  limits  the  precision  of  the 
speed  loop. 

• Decreasing  the  thrust  lag  was  most  effective  for  horizontal 
shears.  For  vertical  shears,  the  aircraft  undergoes  a change 
in  sink  rate  with  no  change  in  airspeed.  Thus,  changing  thrust 
lag  had  a negligible  effect  on  the  Kennedy  shear  hazard. 

• Positive  thrust- to-attitude  coupling  (in  the  sense  of  nose-up 
pitch  with  increasing  power)  has  the  effect  of  reducing  the 
hazard  for  all  three  shears  studied.  Reduction  of  this  coupling 
from  its  nominal  value  (5  ft  thrust  offset)  to  zero  results  in 

an  18  percent  increase  in  Kennedy  shear  hazard  for  the  light  trans- 
port. The  beneficial  results  of  this  coupling  result  from  the 
fact  that  as  the  aircraft  gets  low  and  slow  an  increase  in  power 
to  regain  speed  causes  an  immediate  and  automatic  nose-up  pitch 
which  is  in  the  direction  to  reduce  glide  slope  error.  The 
pilot's  pitch- up  input,  while  probably  modeled  correctly,  was 
not  quick  enough  to  compete. 

B.  TAKEOFF  CONCLUSIONS 

• A suitable  shear  hazard  metric  for  takeoff  is  the  mi  nirmim  rate 
of  climb  resulting  from  the  shear  encounter. 

• For  takeoff  encounters,  the  shears  could  be  ranked  in  the  same 
order  of  increasing  severity  as  for  takeoff,  i.e.,  (i)  boundary 
layer  shear,  (ii)  Logan  shear,  and  (iii)  Kennedy  shear. 
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• The  altitude- dependent  Kennedy  shear  was  capable  of  reducing 
the  rate  of  climb  for  all  vehicles  to  nearly  zero. 

• The  only  practical  operational  technique  examined  which  was 
capable  of  reducing  takeoff  shear  hazard  was  to  reduce  the 
precision  in  tracking  airspeed. 

C.  RECOMMEKDATIDHS 

• Studies  should  be  conducted  to  determine  the  generic  shear 
properties  which  have  a major  impact  on  shear  hazard  for 
various  classes  of  aircraft,  e.g., 


(i)  Investigation  of  range-,  altitude-,  and  time- 
dependency  effects 

(ii)  Spectral  characterization  features. 

• Studies  should  be  conducted  to  develop  the  real  time  use  of 
suitable  wind  shear  hazard  metrics,  given  real  time  wind  shear 
measurements.  This  could  be  accomplished,  in  part,  using  the 
generic  pilot-aircraft  response  estimations  presented  here. 


• Consideration  should  be  given  to  investigation  of  guidance 
modifications  to  minimize  shear  effects  such  as  asymmetric 
glide  slope  sensitivities. 

• Unmanned  simulations  using  pilot  models  should  generally  be 
used  as  a prelude  to  understanding  the  important  effects  to 
be  investigated  in  manned  simulations . 

• Further  work  should  be  done  to  refine  pilot  models  for  i nmanned 
simulations,  especially  where  large  disturbances  are  involved. 
Aspects  which  should  be  addressed  are  trim  functions  (both 
airspeed  and  flight  path)  and  pilot  inattention. 
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APPENDIX  A 
CTOL  PILOT  M)DEL 

An  overall  block  diagram  of  the  CTOL  pilot  model  is  given  in  Fig.  A-1, 
the  key  features  of  which  include: 

• Regulation  of  pitch  attitude  with  control  column, 

• Regulation  of  glide  slope  error  with  a commanded, 
pitch  attitude 

• Regulation  of  airspeed  error  with  throttle,  and 

• Crossfeed  from  pitch  command  to  throttle  for  large 
attitude  commands. 

The  compensation  vised  in  the  inner  pitch  loop,  designated  Y0,  has  the 
general  form  of  a double  lead- lag,  that  is: 

(y+Q  i\r) 

0 °(v*')(w+’) 

The  variables  K(>,  Tj^  , T^,  Tp,  and  Tgff  are  calculated  according  to  rules 
which  will  be  described  later.  This  particular  form  of  pitch  loop  compen- 
sation is  based  on  the  proportional- integral-derivative  (PXD)  model  of 
human  pilot  behavior,  Reference  1 . This  model  states  that  the  output  of 
a human  pilot  in  response  to  an  error  signal  can  be  interpreted  as  the  sum 
of  a rate,  displacement,  and  pseudo-integration  channel  (see  Figure  A-O). 

We  can  write  the  relationship  between  the  two  forms  as 
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Figure  A-1 . Details  of  CTOL  Pilot  Model 
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ERROR 


OUTPUT 


Figure  A-2.  Proportional- Integral- Derivative  (PID)  Model 
for  the  Pilot's  Output  Response  to  a 
Single  Input  Error  Signal 
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The  effects  of  Kr  and  Kj  on  the  system  open  loop  dynamics  are  roughly 
opposite,  in  that  Kr  adds  lead*  to  the  system,  increases  the  phase  for  all 
frequencies,  and  reduces  the  DC  (steady  state)  gain  relative  to  the  high 
frequency  gain,  whereas  Kj  decreases  the  phase  at  all  frequencies,  but 
increases  the  DC  gain  relative  to  the  high  frequency  gain.  Therefore,  to 
achieve  satisfactory  open  loop  characteristics,  we  would  use  either  Kr  or 
Kj  but  not  both. 

For  Kj  = 0 we  have 


where  K0  = ^ and  TL0  = K^ 

For  Kp  “ 0 

Ke(v') 

Y — — ' 

° (T2st1)(Teffs+') 
where  K0  = + K^. 

tl,  * *0  V(kd  * “i)  i t2 

For  either  case,  a crossover  frequency  of  1 .5  rad/sec  with  V?  deg  phase 
margin  was  selected  as  being  representative  of  nominal  (not  maximum  effort) 
pilot  behavior. 

The  procedure  used  to  adjust  the  pitch  loop  compensation  was: 

1.  Set  Tl  = Tl  = T2  = 0,  Teff  = 0.333  sec 


*"Lead"  may  also  be  interpreted  as  pilot  anticipation,  i.e.,  detection 
and  feedback  of  the  rate  of  change  of  a state. 
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2.  Check  the  loop  transfer  function,  Y0  • e/6e.  If  the  phase  margin 
is  less  than  45  deg  at  1 .5  rad/sec,  adjust  Ti^  to  increase  the 
phase  margin  to  45  deg.  Level  Tp  = Tj^  = 0. 

3.  If  the  phase  margin  is  greater  than  45  deg  at  1 .5  rad/sec,  set 
T2  = 10  sec,  = 0,  and  adjust  such  that  the  phase  margin 
is  equal  to  45  deg  at  1 . 5 rad/ sec . 

4.  Set  Ke  such  that  the  crossover  frequency  is  1 .5  rad/sec,  i.e., 

| Y0 -e/se | = 1 at  1.5  rad/sec. 

The  parameter  Teff  remains  at  0.333  regardless  of  the  other  parameters.  In 
this  model  Teff  approximates  the  low  frequency  effects  of  the  pilot's  neuro- 
muscular and  other  transport  delays. 

The  pitch  loop  details,  as  discussed  above,  and  the  procedure  for 
determining  the  parameters  are  common  to  all  three  forms  of  the  pilot  model. 
The  other  key  features  of  the  CTOL  pilot  model  are  (see  Fig.  A-1): 

• A limiter  on  pitch  commands 

• Glide  slope  error  to  pitch  feedback 

• At  long  range,  glide  slope  feedback  gain  is  constant 
in  terms  of  pitch  command  per  displays  (angular) 
deviation 

• At  short  range,  glide  slope  feedback  gain  is  constant 
in  terms  of  pitch  command  per  linear  deviation,  i.e., 
pilot  continually  lowers  his  gain  to  compensate  for 
increasing  display  sensitivity 

• Airspeed  feedback  to  throttle  with  5 kt  threshold 

• Cross feed  from  pitch  command  to  throttle  to  compensate 
for  large  pitch  changes 

• If  the  high  side  0 limit  in  the  glide  slope  loop  is 
hit,  full  thrust  is  applied  and  held. 


A- 5 
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The  values  for  Kd,  Td,  and  were  determined  by  the  following: 

1 . T(j  is  zero  unless  required  to  give  a phase  margin  of  at 
least  30  deg  at  0.3  rad/sec  in  the  glide  slope  loop. 

2.  K,j  is  adjusted  to  provide  a glide  slope  loop  crossover 
frequency  of  0.3  rad/sec  at  short  range. 

3-  Ku  is  adjusted  to  provide  an  airspeed  loop  crossover 
frequency  of  0.15  rad/sec  (neglecting  the  feedback 
threshold) . 

The  airspeed  feedback  and  flight  path  command- to- throttle  crossfeed 
thresholds  were  determined  by  review  of  previous  piloted  simulation  time 
histories.  The  values  were  ±5  kt  airspeed  and  +5  deg  path  angle,  respec- 
tively. 

The  crossfeed  gain,  Kqt,  is  an  acceleration  command  to  the  throttle  to 
compensate  for  large  pitch  attitude  excursions,  and  thus  has  the  value  of 
g,  32.2  ft/sec2/ rad. 

The  speed  loop  gain,  K§T,  is  the  inverse  of  the  x-acceleration  deriva- 
tive with  respect  to  throttle  (i.e.,  KgT  = l/XgT).  It  is  simply  a convenient 
way  of  scaling  the  acceleration  command  to  a throttle  command. 

Tables  A-1  and  A-2  present  summaries  of  pilot  model  parameters  for  the 
three  CTOL  aircraft  used  in  this  investigation  as  well  as  for  the  STOL 
commuter  aircraft.  The  STOL  pilot  model  is  discussed  in  the  text.  Section 
II .A. 2,  and  the  takeoff  pilot  model,  in  Section  II. A. 3- 
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AIRCRAFT  MDDEL  DESCRIPTIONS 


The  following  tables  describe  the  aircraft  models  in  terms  of  fixed 
operating  point  linearized  equations  of  motion.  Stability  derivatives 
and  transfer  function  factors  are  given  for  the  landing  cases  (Tables 
B-1  through  B-4),  and  only  stability  derivatives  are  given  for  the  takeoff 
cases  (Tables  B-5  through  B-8). 
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~£ 

.61  O 

8 . 07  >C  - 

,16025.219  -.03 

50  5 .216  ] < 1 . 007) 

Mb 

THE- 

-DTH  > 

• 

1212C 

.8345.3 

42  .285  I . 1885  1 

.01415) 

iNi)  ' 

.<  DP- 

-PTH  ) 

■5 

he 

.61  < . 

1 393  ’C  . 

55752.56  1.4255 

2.13]  <2.39 

(Nl*)  "•*=■ 

-PTH' 

“C. 

.61  < . 

010  30 ' ' 
-.0238.' 

.13)8)1  .553  52.57 

1 .42352.14] 

(/Si,)  M<  HD- 

-DTH  • 

I 

.46.2' 

. 1507  )t 

.662  5 3.36  2.22  5 

2.52]  <2.49) 

v 

Numbers  enclosed  In  parenthesis  are  first  order  factors,  (s+a)  » (a). 
Quadratic  factors  are  enclosed  in  brackets,  ls2+2£  <ue+cuP|  . U;  ujt.  Other 
numbers  are  the  transfer  function  gains • 
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TABLE  B-1  (Continued) 


r . ■ 


PE  PTH  rOUEL IMG 
(N^v)  M'U-PE  U-DTHi 
(^fc  i^)N'  U-PE/'THE-PTH) 
C^hetr)  h-U-I'E . PD-PTH) 
N<U-PE/fiZ-PTH> 
N<U-DE'HP-PTH> 
M(W-DE-THE-PTH> 
rKw-DE/rr-HTH) 
M'W-OE'flZ-PTH) 
M<M-PE'HP-PTH> 
M<  THE-PE-  PP-PTH) 
M< THE-PE  0Z-PTH) 
M< THE-PE -HP-PTH) 
H(PP-I'E-  0Z-PTH) 
N*  PP-BE<HP-DTH) 
fKRr-BE'HP-DTH) 


MIJMEPHTOB : : 

157. 2<-.  0214X14. 56)  <-48. 9) 

9.55<  .461  > <4.40 
-157.2<2.63X-2 .67)  < It  04 . > 

157 .2<  .00662X2. 65K-2. 69)  <-7.4! 
-156.9<2.63><-2.67>  <1102.) 

-2. OKI  .251  ) <-2.51 
-462. <1. 251 > < -6  03  . > 

482. < .00703X1  .251  > <4.23> 

8 . 22  < 1 .6 1 3 X —45 .6 1 <-604. 

-2.01(1  .251)  <-2.5K 
2.0! < 0.0)' 1 .251 > 2.51 ) 

-1 . 506< 1 .513)  < -2 .28 > 

-3.38(1.251'  <-4.23> 

-8.22(2  .63> (-2 .67 1 <57.6 
8.22<-.0696X-2.57X2.6l  > 3.84> 


■77*» 

» 

r.\ 

m 

c 


I .nr.  NIJMEEHTDE'S  : 
(N  uj)  M(U-I.M5> 
M(W-IJG) 
nt  the -mi;  > 
n< pp-ug> 

N(HZ-I.'G> 


.0337(  .871>t  .254 : .685  .17421.663  3 <.01536> 

.257 ( .876E-7X  . 001431  X .480)  < .759E-11> 
■464E-4(0. OX -10.23)  <-.000474 : 

-.257<  .186  3E-8)l  .3295.666  .2191.6£93  <.?04€-tl> 


2151.6273 

M(HB-UG)  - .258<  .00703)t  .3285 .665  .2181.6263  <-.000604> 


357(-.169tE-6 X .00707)1 .3251 .66' 
< .330E-12 > 


Mi;  MUMEEHTOPS  t 

H<  ij-mi;  > 

Md.l-i.IG ) 


- . 0479<  . 1863E-8  ’1  -.  0844  1 .694  -.07555.8913 

<- . 1302E-9) 

.03271  .07805 .1501  .01 1831 .1496  3 

1.67354.57  3.0713.383  <.01536 


M<  THE— WS  • -.001  366<  .373E-8X  .021  3)  <-.6  56  > <-.44  0E-12> 


n- 

I'D 

-W6> 

- .0  27< 

.834)<  15 . 

38  '1  - 

-.015  36 

f .I 

mZ 

-MS  > 

.529' - 

.559E-8' < 

.802) 

1 - 

.03931.1900  - 

m* 

HI'- 

* W13  > 

- . 0 ' 26< 

.84 : i5. 

30)1  - 

-.01534 

B-3 
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TABI£  B-1  (Continued) 


PE  IJG  COUPLING  CUNEFflTORS : 

N<U-PE  l.i- UG > ,£4£<  t . 113X13.42)  <3.61> 

mU-PE  THE-UG > . 01469<  .847)  '.01244') 

mu- PE  PP-IJG  - . £4£C  - .01131  51  .613  - . 01823  5 1 .6  t 2 3 <- .629) 

NCU-PE'EiZ-IJG  > .242<-.0  847>t-.000555»l .580  -.000878  51 .580  3 

<-.081 0> 

N<U-DE -'HP-UG  > - .2481  -.  01 131  51  .61  3 - . 0 1823  51  .618  3 <-.628> 

M. u- PE  THE -MG  ' .11181 0.0)  .1116) 

mW-DE<DD-UG>  £6 .8<  .259E-9 > < .693E-8) 

N<W-PE  8Z-U6)  -£6.8<  .444E-10X  .00703)  <-.836E-U> 

MO.I-PE  HP-IJG  > . 01266 < .00703)  <£130 . ) <.1891> 

M*  THE-PE/"  PD-UG  1 . 1 1 16<  — . 749E- 10)  <-.836E-ll> 

M< THE-PE  8Z-UG  > - . 1 1 16< 0 . 0 J <- .749E-1 0)  <.836E-11> 

N<  THE-PE-- HP-UG)  .1 1£2< .00580)  <. 000851) 

M<  PD-BE/'ftZ-MG)  . 1881  < .341E-9  > <.64£E-10> 

M<PP-PE  HP-UG)  -.01£88t-. 0113151. 813  - .01823 5 1 .81£ 3 <-.03£9> 

h < HZ- PE  HP-IJG  ) . 0 1 £ 68 < . 0 07 03 X 3 .49 X -3  .53 ) -.001  098 > 

PE  'U6  COUPLING  NUMER8TDPS s 

MXI-PE  I.I-U6)  -.355<-l  .350X13.71)  <6.57> 

rUIJ-PE- THE-WG)  - .0£04<  .1  382E-8)  <-.282E-10> 

N<u-PE'PP-UG  ■ . : 55<  3 .65X-5 . 06)  <-6.57' 

miJ-PE'NZ-WGi  - . 355  < .18  39E-8X  3 .66 X -5 . 08)  <-  .£8£E-9'> 

N>U- PE  - HD-146  i .354(3.65  '<-5.06)  <-6.58> 

M <i.i- PE  - THE-i.n?  > .0£40<  .0608X8.52)  < .01244) 

mU-PE'PD-l-IG  ■ 5. 76<  .0608X8. 5£>  <2.98' 

mi.i-PE ^82-1.18  I -5.76<  . 00703 X .08  08  X8. 52)  < - . 0210) 

N<  1.1- PE  HP-1.16  ) 5.74c  .0<Gi.X8.50>  " ?.3£  ■ 

M<  THE-PE^DP-WG ) .0240<  .0608 X8 .52)  <.01244 

fl'THF-PE  8Z-MG.  - . 0240*  0 . 0 ' < . 0608  '<8.52)  <-.01244) 
m THE-PE  HP-i.»J  > .0240*  .061 1 '<8.47)  <.01243' 

h*  PP-PE-HZ-i-'G  .0405<  .0608X8 .52)  .0210 

m I P-PE  <HP-l.'G  > . 0 1 656  < 3 .65  X -5.06)  <-.344' 

ft'HZ-PE  HP-i.ii"..  - .018*8  < 1 . 3!9>f  - . 312!  .918  -.2665.8723  <-.0209) 
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TABLE  B-l  (Concluded) 


DTH-IJG  COUPLING 
H'U-PTHxW-IJG  ■ 
mU-BTH'THE-IJG ' 
NCU-PTH-  PP-IJ6) 
rKU-BTH-'PC-UO) 
IXU-DTH  HP-UG) 
N<  M-DTH'  THE-UG  ’> 
mW-DTH^PP-l.ii?  ■ 
N<U-PTH''8Z-UG> 
N<.U-BTH- HP-UG) 
NCTHE-DTH'PP-UG  ' 
N(THE-BTHrflZ-U6> 
Hr THE-PTH^HD-UG> 
M< DD-I'TH-F,Z-UG) 
N<BB— BTH-’HB— UG  > 
M«  RZ-PTH'HB-IJG) 


liunEBftTOBt : 

9.71«-.£91  X .606)  < - 1 . 007-  > 

-.00307<4.6I  <-.0141* 

-9. 7U  .3435.646  .£££5.607]  <-£.39> 

9. 7K. 01008>t  .3375.644  .£  1 7 8 .GO*  ] <.0£.‘8> 

-9. 7U  .3435.648  .£££5.607]  <-£.38> 

- . 0 3 1 £ ' 0 . 0 > <-.03!£'> 

-7.4-'  ,£40E-9  > <-.180£E-8> 

7.49<-.£67E-9X  .00703)  <.- . 1 408E- 1 0- 
.£99< . 0071 3 ' < -£4 .7)  <-.09£7) 

-.031£<0.0)  <-.031£> 

.0312(0. 0X0.0)  < .03t£' 

-.0314<  .0£36>  <-. 000741) 

- . 09£7<-  .£67E-9  ) . 1408E-1 0 

-.£99[  .:4.?5.646  .£££5.16.07]  <-.l£49> 

.£99<  .00706)1  .£845.7169  .£185.737]  <.001£48) 


PTH-  UU3  COUPLING 
N<  U-PTH--W-WG) 
ri<IJ-PTH  THE-WG) 
N<  U-PTH  PB-WG> 
N<U-BTHrfiZ-MG-. 
U'lJ-BTN/HP-WG) 
N(W-PTH  THE-HG) 
N(W-PTHrPP-UG) 
N<W-PTHr«Z-UG> 
N'I.i-MH  HP-WG) 
N< THE-PTH- rp-HG) 

n> the-pth-hz-ug> 
rt<  the-pth-'Mp-i.'G  > 
N<  PP-riTH-HZ-'.'G  ' 
ruPD-DTH  HP-W6) 
fKHZ-OTH-HD-UG  > 


MUMEPHTOR'C  : 

• 719<-.1£37>C  .63354.80  ; .045  3.7£]  -£.04> 

- . G£99<  - .709E-8X-1  .068)  <-.££9E-9> 

- .719< - . 1937X  .968 ’ < 19  .££)  <£.04> 

.719*  0 .0)  < - . 1 97 1 ) (.  .978) < 19  .£t ) -£.09> 

- .714<-  . 1937X .968)' 19.££'  <£.04> 

- . 000 389  < .£4  ?X 149 .6)  <-. 01415  > 

- . 0933  < .£43X149.6)  <-3.39- 

. 09 ?3<  .00714  X .£43X149.6)  < .0238) 

.0  :74<  .£66 X -17 .££X£0.4)  <-3.90' 

- .000389<  .£4  3X149.6)  < - .01419) 

.000389' 0.0X .£43 X 149.6)  .01419) 

-.0019*4*  .£9£X£8.7)  -.01413) 

-.000696*  .£43X149.6)  <-.0£38> 

- .0374*  - . t9  37  .968  X 19  .££)  < .106-7) 

. 0 ";74<  .893) < 19  .££>C  - .1902  5 .££1  -.04£1  5 .£17  ] 

7 . Or  : 8 
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TABLE  B-2 

STABILITY  DERIVATIVES  AND  TRANSFER  FUNCTION  FACTORS 
(LIGHT  TRANSPORT,  IANDING) 


GEOMETRY: 


VT 

GAMMA 

219.8 

-2 .994 

A 

PHD 

MACH 

XIO 

ZJ 

1118.7 

.002377 

• 1 9667 

3.390 

4 .999 

£ 

C 

WEIGHT 

IY 

ALT  I TIJDE 

9:- 0.0 

1 1 .200 

97000. 

877100. 

0.0 

NON-D 1MENS ICNAL 

DERIVATIVES: 

CL 

CLA 

CLAD 

CLM 

1 .8799 

7 . 097 

-8.000 

0.0 

CMA 

CMAP 

cue 

CMM 

-1 .4724 

-2.400 

-24.97 

0.0 

CP 

CPA 

C DM 

TM 

TPTH 

.2708 

. 7899 

0.0 

-3720. 

26840 . 

CPPE 

CL  PE 

C MPE 

0.0 

.4783 

-1 .6700 

DIMEMSICNAL  DERIVATIVES! 

XU 

XU  * 

XW 

TO 

-.04  347 

- .04471 

.07708 

-.0010676 

ZU 

ZU  * 

ZWP 

ZU 

-.2910 

-.2909 

.017671 

-.6362 

MU 

MU  • 

Mwp 

MW 

-.0004773  ■ 

- .0004934 

-.0002004 

-.004620 

MAP 

MA 

Mp.j 

— . 04  399 

- 1 . 0783 

- .4774 

::pe 

ZPE 

MPE 

0.0 

-8.716 

-1 .1879 

XDTH 

ZDTH 

MDTH 

9.074 

-.7  377 

. 17337 

DENOMINATOR : t. 09675. 1614  . 01593  5 . 1606  ] 

I. 49451. 177  .781  51  .023]  <.0361> 


DE  NUMERATORS  : 


N1 U-PE  ) 

Hi  w -de) 

N« THE-DE) 
N<  DD-DE) 
N< AZ-PE  ■ 


- .6641  I . 273)  1 — 26  .6  ) .23.2 

-A. 87130. 3 >t  .12405.204  .02731. £03]  <-lt.£4> 

-1 .184< .G877X .789)  <-.0790> 

8.87<  .0I093X- ■ .71X4.77)  <-l  .716': 

— 8 .87<-:3 .74  X4 . 79  >C  .381  1.0270  .009735.0231  ] 

.099: 


N'HD-DE)  8.66<  .01874X-3.72X4.78)  <-2.92> 


DTH  NUMERATORS  I 
N'U-DTH  ■ 
N< U-PTH) 
N(THE-PTM) 
N< PD-DTH' 
N* hZ-PTH) 


f"  HD-PTH  > 


9.07<-,£4t)I  .77751  .208  .8935  .988]  <-3.18 

-.7475-77.2'!-. 03145. £17  -.008825.217]  t.473> 

.1737t  .9®4  5 .344  .1:425.0387]  <.01811) 

.747<  .1078>C  .32858.71  2.1258.18)  <2.70> 

-.7471  .OI282X  .10t  3H  .32358.77  2.1258.20] 

- . o ; 07 

.07271  .117' ,c  .712517.89  1 2 .79 5 12 .43  ) <2.G7> 
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TABLE  B-2  (Continued) 


DE  DTH  CCUFLIM6  NUHEPhTORS s 

NUJ-BE- U-BTH)  SO. V- .01890X30.4)  <-46.2> 

fUU-BE^THE-BTH)  10.75<.598>  <6.43> 

H<U-BE  BB-BTH)  -80 .5< -3  .77X4 .80)  <1457. > 

M<  •J-BE''FC-BTH>  80 ,5<  . O0732X-3 .80X4 .82)  <-10.81> 

N<U-D£- HD-BTM)  -80 .4<-3 .77X4 .80)  <1455. > 
fKM-BE  THE-BTH)  -2 . 01  C l .61 1 ' <-3.24> 
rUH-BE'-DD-BTH)  -441.  <1. 611  > -711.) 

M<W-DE-RZ-IiTH)  441  . < . 00788  X 1 .8 1 1 ) <5.44 
N<W-PEHP-BTH)  4 .2! <8 .2 1 X -76 .7)  -712. 

N<  THE-BE  DD-BTH)  -2. OKI. 811)  <-3.24 
N<TH£-BE'«Z-BTH>  2. 01 <0.0X1 .811)  <3.24' 

M< THE-BE -HB-BTM)  -1.448(2.00)  -2.90> 

H<BB-I'E'8Z-BTH>  -3.38(1.611)  -5.44> 

N<BB-BE^HB-BTH)  -4.21 < -3 .77 X 4.80)  <76.1 
M«flZ-BE'HD-BTH)  4 .21  <- . 0664X- 3 .67X4 .76)  <4. 87) 

IJ6  MIJHER8TCRS : 

M<  1.1-1.16)  .0445<.989>C  .3721.905  .337 5 .340  J < . 0361  > 

n<l.l-IJG)  .296t- .596E-7X  .00345X  .826)  <.67tE-10> 

M'  TMF-U6)  .0004  34<  0 .OX-2  .55)  <-.001106 

WPB-IJG)  -.296<-.1863E-8'I  .2805.908  .254S.872J 

- .296E- 1 0> 

M'RZ-IJG)  .296> -.987E-7X  ,00769)C  .2775.904  .2505.8691 

< — • 926E- 1 2 > 

N'HB-IJG'  -.293<  .00765)t  .2805.909  .255  5 .872  1 <- . 001863) 

1.18  MUMERHTOF  3 • 

H<U-Hi3)  -.0785' 0.0 'C  .0599  51  .209  .072451  .2071  <-.1400E-8) 

M« M-M6 ' -.01 799 ( 4 . 62  X - 1 6 . 77 ) t . 1 447  5 . 1 6 09  . 0233 5 . 1 592  1 

< . 0 36 1 > 

N<THE-MG>  - .001879«  .67tE-7)<  .0421  X-t  .891  > ..1852E-11> 

fi.  DB-i.tr-  . ,01799<  .986  .<-36.0  >t  .00703  5.238  .001670  5.2381 

-.0361 

N<«Z-'.«5>  .60<  0.0X1  .025  C -.010045  .235  - . 00236  « .235  1 

. 370E-I1 

H*  MB-iG  . .01  7°5  < .999).  -35.8  >[- .01 1 041 .237  - .00261 5 .237] 

- . 0 > 0 
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TABLE  B-2  (Continued) 


i 

PE  UO  I tXPL  IH6  HURCRRTOB'  : 

M<  l.i-BC'M-US'  . ,97ft  ■ 29.4  > <11.£4> 

M*  L'-PE- THE-UG>  .0"527<  1 . 120)  <.0'590> 

N'l'-PE  PP-IJG)  -.39V-1  .632X2.66'  • 1.716) 

M'I.i-pe  RZ-U6 > .29ft  .OM4X-!  .717X2. 69  > (-.0993) 

N(U-be/,hb-ug>  .6;:2><g.66>  <».ri3> 

H<  U-PE  - THE-'JG  ■ .347(0.0/  (.347> 

MU-PE- BB-'.'G  ■ 76.£(-.286E-9>  <-.£16E-7> 

M<  IJ-PE  RZ-UG)  -76. 2( .1206E-9X .00766)  C-.704E-10) 

N(W-BE'HP-U6>  .0£06(  .00766X3720.)  < .^87> 

N<  THE-PE- PD-IJ6)  . . 47<  . 1 068E- 1 0 > .370E-11) 

N(THE-PE'RZ-UG>  -.347(0. OX  .1P68E-10)  (-.370E-11) 

N(THE-PE-HD-'J6>  .349(. 00863)  (.  00308 > 

THEE- PE' RZ-UG)  .683<  . 1206E-9)  (.704E-10' 

M<  PP-PE-HB-UG)  -.0£06<-l  .632X2. 66)  < .0686) 

M<RI-BE-HP-IJG)  ,0£06<  . 00766 )( -3  .24X6 .£6)  (-.00518) 

PE  -'MG  COUP L IMG  HUHERRTCPS: 

MdJ-PE/i'i-WG)  -.684(-l  .071X31  .6)  <£3.£> 

M(U-PE/'THE-l-'G  ' -.0'-*16(-.16I7E-8)  (.1482E-9) 

h<IJ-DE'PP-U6  ■ .684<-5  .68X6.41  > -£3.£> 

nnj-PE  RZ-UG)  -,684(  .691E-6X-S. 30X6.43)  (-.237E-9) 

MCJ-PE- HD-WG)  .68  3' -*•.28x6.41)  <-£3.1  > 

tt<  lii-PE-  THE-WG)  -.0046?(  . 08  1 4X -1*56 .7)  (.0590) 

N'U-DE  EP-MG)  -l  .016<  .0814X-156.7)  (12.96) 
r • u-DE  RZ-U6)  1 .016<  .00760X  . 081 4) (-156 .7)  (-.0993) 

N<  I.I-PE/HP-MG)  -1  ,051<  .0896X-1  *0.4)  <14.15 


tt<  THE-PE.  BB-WG  > 
M< THE-PE  RZ-WG) 
f«  THE-PE  HP-WG) 
M< PP-PE/RZ-HG) 
M' PP-DE'HB-WG) 
M<  RZ-PE' HB-HG ■ 
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dthug  i:n"PLtns 
M<  U-PTH<U-UG  > 
M*  U-DTH - THE-U6  > 
fKU-PTR'DD-UO  • 
N<U-DTH-BZ-U6> 
H(  U-DTH  'HD-UG  > 

N<U-DTH-  the-ug) 
fKM-DTH- DD-UG) 
N( U-DTH- 8Z-UG  ■ 
M<  U-  PTH<  HP-UG  > 
M<  THE-DTH- DD-U6 > 
N<  THE-DTH- 8Z-UG) 
M< THE-DTH- HD-UG > 
fl<  DP-DTH- HZ-UG  > 
M<  PP-PTH  HD-U6) 
H«  HZ- PTH  HD-UG) 

DTH  -Ui;  C DUEL  IMG 
H< U-DTH  u-WG> 
MlU-DTH  THE-WG • 
M<  U-PTH'DD— HG) 
U-DTH- HZ-UG  > 
M<  U-DTH-- HD-UG > 
I1(U-DTH'THE-W6  > 
M<  U-DTH  PP-UG 
M<  M-PTH'RZ-MG • 
M<  U-PTH<HP-UG> 
N(THE-PTH<PP-UG' 
H< THE-DTH  HZ-l-G) 
ru  THE-DTH  HD-UG) 
H'  PI -PTH  HZ-UG’ 
H> PP-PTH  HD-UG- 
fi<HZ-DTH  HD-UG  > 


B2T  AVAIL JM 


TABLE  B-2  (Concluded) 


MUHEBBTDH’ t 

2.711-.617X  .881)  (-1.473) 

-.00289(6. £6)  <- . 0181 1 > 

-2.71C  .2398  .*6t  .249 ! .928  3 <-£.30> 

£. 71  < .01£2£>C  .£33  8.939  .2438.927]  (.0303) 

-2.71I  .£39!  .961  . £498.928]  <-£.30> 

- . 0437 (0.0)  (-.0437) 

-10.031-.il 34E-G > <.1138E-7> 

1 0 .031-  .771E-9) 1 .00766  • 1-.393E-10 
.1416« .00771 >1-70.3)  <-.0769 
-.04371-. 649E-9)  <.296E-10> 

.0437(0.0  X-.649E-9)  <-.£96E-10> 

— . 0438  ( . 0£07  ' <-.000946) 

- . 0768<  .386E-9 ) 1-.296E-10) 

-.1416r  .£398.961  .£498.928]  (-.1308) 

.1416*  ,00768>C  .£0381  .210  .£4381  .183]  <.001391 


NUHERHTOPI  : 

- . 1632* - .£04 ) <3 .34) <- 17 .£3)  -3.13) 

-. O1703<0. 0)(-£.-0)  <.0444 
.16 3£<  — .4 1 1 ) ( 1 .3£7'<-33.7>  3 . 18 

-.16321-.286E-4X-.41?>*  1 .337)  ( -33 .7)  <-.474E-9 
.1630"  - .41 1 )<  1 .327U-33.7)  < 3.18> 

.00379*  . 1 836 ) < -23 . 7)  (-.01611 
.632*  .1838)*  -£5. 7)  <-!  .98) 

-.6321 .00767)1  .1637X-23.7)  < .0303' 

- . 0 063  3 ( .£01 >1-30. 7) 1-78. 9)  (-4.14) 

.00379* .1838)f-£3.7>  (-.01611 > 

-.00' 7910.0)'  .1638) (-£3. 7)  < .01811 > 

.002891 . 1906)1-  32 .6)  -.01809) 

.006371 . 1 638 ) ( -£3 .7)  (-.0303) 

.0063  <- .41 1 ><  l . 32?><-33.7>  .1663) 

- .0083  Ml . 0 ! 0 >1-  33 .7 't  - .1636  1.311  -.  0376  8 .306  j 

* .0304- 
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TABLE  B-3 

STABILITY  DERIVATIVES  AND  TRANSFER  FUNCTION  FACTORS 
(BUSINESS  JET,  IANDING ) 


GELMFTRYi 


VT 

. .hum* 

£11.1 

- .000 

R 

PRO 

**CH 

XIO 

ZJ 

11  If.  -c 

. 002377 

.1291  0 

6.071 

-1  .4979 

c 

c 

ME IbHT 

1Y 

*LTITUDE 

c*:i ... 

7.040 

12300. 

18701  . 

0.0 

MOM-rirEh'IOMRL  DERIVftT  IVE 

C ; 

CL 

CL* 

i L*P 

CLM 

I .0001 

4 .332 

1 .7600 

0.0 

CM* 

CM*P 

c*0 

CM* 

- . ^<*69 

-4.200 

-12.000 

0.0 

cri 

CD* 

CD* 

TH 

TPTH 

.17104 

.6163 

0.0 

0.0 

39 . 00 

CI‘PE 

' LPE 

C *PE 

0.0 

.4406 

- 1 . 0666 

PIMEMCICMhl  PERI V* FIVES* 

::u 

XU  ♦ 

XW 

TIJ 

- .091 £2 

-.oci£2 

. 03742 

0.0 

cu 

ru  ♦ 

7WP 

ZW 

- .£-994 

- .£994 

.004393 

-.7070 

MU 

MU  ♦ 

MI.IP 

MM 

.0011081 

.0011081 

-.017321 

-.02183 

flRO 

MR 

Mi) 

- . 3899 

-4.608 

- .2248 

XDE 

ZT'E 

MPE 

0 . 0 

-13. *24 

-3.032 

I'TH 

ZDTH 

mpth 

.1503? 

- . 0£  1 3£ 

- . 004726 

I>EMGMIhRTO»r  C.  12 15 5 

.£  '>■ 

.21?  .0239 

5 .211  3C  .43 

752.28  . 

I«E  MUMF PRIOR 

o - 

mu- 

PE)  -,796<1 

.03351- 126. 

15  1 03 •?> 

m i.i- 

PE  ■ -1  .6  6 < 

77.2  51  .136' 

5 .214  .02935.2123 

n< THE  - 

DE  ' -3.01  < . 

08205*  .616) 

<-.253/ 

ru  ri- 

DE>  13.86< 

.00628)1-6 . 

49 ' < 7 . 3 0 > 

s-4.24) 
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hZ- PE  ) -1  3 .66  < -►  .'4  )<  7 .345C  .297  5 .0230  .00742  5.0239  3 

.426 

M<  HP-DE  i t:.  4<  .01422M-6.32  >17.32)  <-9.66> 


12'. H 2.07DE6  24.767.  * C=LF* 

PTH  HUMERfiTORS* 

M»  l.i— Ii  TH  ■ . 13P4<  .1 14'-  C .395*2.  28  .93712.123  <.0970 

fu  M-  OTH  > - .0212<  4'-*  .2  >C  . 1476: 1 .209  .03095.2073  - . 0462  ) 

M'THf  DTH>  - .00464'  - .222X  .832>  <.000667) 
re  PP-  I'TH  I . 02121- .4'7><  - 3.  33  ><  7.06)  <.229> 

fl'*Z-PTH>  -.0212'  .00641  ■ < - .446  '<-3.39X7.09’'  <-.0014615 
ri'Ml-t'TM'  .00  <- .400  " -4 .64"  9.04''  <.224> 

>-3  B-10 


TABLE  B-3  (Continued) 


PE  /DTH  COUP Lirti?  MMMERHTDRS : 


IKIJ-PEU-DTH'  2. 08c-. 000363X77. 2)  <"-.0584) 


fi<U-PE  THE-PTH' 
MOJ-PE-'PP-PTH) 
M<  U-PE ' HZ-PTH  > 
mU-BE-HD-DTH) 
h' 1.1- PE  THE-PTH) 
MCW-PE/ DP- PTH > 
M-  W-DE-'HZ-PTH) 
McM-PE'HP-PTH) 


.753<.t44>  < .485> 

-2 .08<-6 .54) <7 .“"to  <102.4/ 

2.08<  .00788X-8.58X7.54)  -.817) 

-2 . 08  c -6 .54) <7 .50)  <102.2> 
-.0413<5.5?)  <-.230> 

-8.72(5.57)  <_  -48.8) 

8.72<  .0078v- X5. 57)  (.388) 

.lO-Xc  19.83X-22.4)  <-48. 5) 


M- THE-PE/DD-PTH) 
ti • THE-IE  8Z-DTM) 
rKTHE-DE/HD-PTMT 
N- PD-PE'HZ-PTM> 
N< PD-DE  HP-PTH> 
McflZ-PE  HPc-PTH) 


-.0413(5.57)  <-.230) 
.0413(0.0X5.57)  <.230> 
-.001832(111  .5)  <-.204> 
-.0888(5.57)  < - . 388 > 
-.1081C-8.54X7.50>  <5.36> 

.1081 < -.0844) (-8. 51 X7.53)  < .344) 


IJG  MUMFPflTDlX  : 

TiC  U— 1.1(3  * 
Mi  M-|JG> 
ft(  THE— IJG  :• 
Mi  PP-IJG) 
Ml  HZ -Mi-  ■ 


.0512(  .818)t  .31512.24  .708:2.133  <.236> 

.288*  - .871E-7) (-  . 0350X  .1784>  < .l70;E-8> 

- .0018  301  0.0X4.47)  <-.00728 

-.248<  .261E-7H  .28852.27  .84852. 183  < .827E-9) 

.288c  .81 1E-GX  .00798)1  .28452.27  .84552.181 

< . 181 IE-1 1 > 


MCHP-IJG)  -.300(  .00798)1  .288  52.27  .65052.18)  <-.01238 


12-XET  2. 0'PEG  24.76>.  ft  C=LEO 
MG  mjMEE8TD(-t : 

M'U-M6>  -.0572c  .233E-7)!  .056653.24  .183453.23) 

<- .475E-8) 


M'U-MG ' 


M« THE -MG  i 


ro  pr— mg) 


.004371 .1236  5.212  .0282  5.210) 

1.513534.7  17.78529.7]  < .236) 

-.00284c  .536E-8X . 0438 ) c -7 . 07)  < .724E-11) 

- .004  371  1 .929  >c  160.9 >1-  . 01288 5 .417  -.005295.417] 

- .2  36) 


M>  8Z-U6  | .712c -.  1863E-~  •'  1 .934> 

1-. 01775. 41?  -.015555.412)  <-.£90E-10) 

M>  HP- MG  - .004371  1 .9*  8 ■«  160 .2)1 -.  0385  5 .415  -.015985 .415  ) 
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TABLE  B-3  (Continued) 


PE  IJG  'DUELING 
M<  l.l- DE  ■ 1.1-1  nV> 
MXJ-PE  TME-HG  .■ 
N<U-DE  PD-UG) 
N*  IJ-DE'HZ-U'3  • 


fiUMf.PPTriPS  : 

.710<  .907  <7*. 3)  <49. £> 
.£5E< .988 ' .£53> 

- .71  Of-*  .01  ><£  ;9?>  <4  .£4 
.710i  .090£X -£.18)<3.G5> 


<-.4£E> 


M<0- PE  HD-UG 
N<W-PE  IHE-l'G 
M«  U-DF  DD-UG > 
ncg-PE  82-UG) 


-.708C-e.OncE.f7>  <4.£3> 
l .Ctv  0 . 0>  < 1 .515  ' 

'3£0.<  . 1449E— 9)  < .484E-7) 

-3£  0 . < . £89 E-9 ' C . 0 0798  > < - .7  39E-9 > 


N<  U-PE  HP -IJG  ) 


M' THE -PE  DD-UG  > 
MUHE-PE  HZ-UG  > 
M'THE-DE  HP-UG > 
M<  PP-PE  -HZ-UG) 
MCPD-PE  HP-UG 1 
N<HZ-PE  HP-UG > 


.037EC  . 00798x8670. > <£.57> 

1 .515C-.GG1E-10>  <-.1304E-9> 

-1  . UJ 1 5 < 0.0)0  .8E1E-1  0>  < . 1 304E-9  > 

1 .5£7< .0088-7)  < .01384 
£ ,55<  ,£89E-9)  < .739E-9> 

- .0'7£<-£.01)<£.97)  <.£££> 

.037£c . 00738 >< -8 .£1X9.18'  -.0££3> 


PE  MG  CGUPL  IMG  NIJMEMT0R8  t 

M'U-DE  IJ-WG  ) -.798.C-1  .701X74.8)  <103.7) 

M'  IJ-PE'  THE-MG  > -.£84<-.16  £E-8>  .4E4E-9) 

N*  l.i-PE-  PD-WG  > .798  ' 1 1 . 16)* -t  1 .67)  -10:.7> 


MC  iJ-PE  HZ-Mi?'  - .7*0  .M1E-6X1 1 .19X-11  .70)  .£3£E-9> 


M< U-PE  HD-WG ' 
M<  1.1- PE  ^ THE -MG  ) 
H<U-PE  DP-MG  > 
H*  l.l-  PE -HZ-MG 
Mi  1.1- PF  HP-HG) 
M‘  THE -PE  PP-i.iGJ 
M<  THE-PE  HZ -MG  > 
IK  THE -PE  HI'-MG  ■ 


.795C  1 1 . 1C  ><  - 1 1 .G7>  <-103. 
.0841 .0781X55.4)  <.£53> 
1£.  4<  .0781)155.4)  <53.4 
- 1 £ . . 4<  . 00787 X . 0781  > < 55 ,4> 
1£.£8<  .0885)155.3.)  \58.7> 
.05841.0781X55.4'  .£53> 

-.0e84<0.0 '<  .0781  X55.4)  <- 
.0584*  .0785X55.1  • .£53 


5; 


<-.4£G> 


.£53: 


M> PP-PE  HZ -MG ) . 0985< . 0781  • 55 .4)  < .4£4> 

fl.  PP-PE  HI* -MG  • . 0417C  1 1 . 18  •< -1 1 .8  7)  <-5.4 

li<8.  -PE/-HP-I.G  - . 041  7*  £ .97  H - . 301  it  .858  -.55851  .7701  -.4£8 
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TABUS  B-3  (Concluded) 


FTH/UG  COUPLING 
M<U-PTHU-IJ6  • 
M<U-DTH'THE-UG) 
fi-  U-I'TH  DD-UG ) 
N<U-PTH-ftZ-U6> 
h<  u-pth-'Hp-og  ) 
M<W-DTH  THE- -US' 
M<  M-DTH  DD-UG) 
h<  W-DTH-  hZ-UG  > 
M<  M-PTH  HP-U6) 
N( THE-PTH  rD-UG) 
f1<  THE-  DTH-'HZ-UG  ) 
N<  THE-PTH  HP-UG) 
M(  PD-DTWiflZ-UG  > 
Ml  DP- PTH'HD-Ub ) 
M<6Z-DTH'HP-U6> 

DTH  WG  C0UPLIM6 

MU.l-PTH'-H-M?) 
ft'  M-PTH'- THE-WG ) 
N'  U-PTH'DP-MG) 
M<M-DTM)HZ-WG) 

f1<M-DTHHP-W6» 
H'  U-DTH'THE-WG) 
H> M-PTH'PI-MG) 
M<  W- PTH  HZ-U6) 
M<  W-PTHiHP-MG  ' 
M' THE-PTH'  DP-IJ6  > 
Ml THE-DTH'PZ-MG ' 
Ml THE-PTH  HP-WG» 
H<  PP-PTH'HZ-MG i 
M<  DP-DTH'HD-MG • 
H<  HZ-PTH-HD-MG  • 


MtJMERBTnP?  : 

.04591  .63' 51  .003  .6335.776)  1.0468) 

-.4«E-5(t73.8!  <-.000667 

- .04*91  .£69  58.83  .'■47’ 58 . 14  ) <-.2£9> 

.04c9< .00636)1  .£8882. £3  .64358.14)  <.001461) 

-.045,-*t  .889  58.83  .64758.14)  <-.829> 

.001438<-.63£E-9)  <-.906E-l£ 

• 30£<  .38  3E-9)  <.1159E-9> 

- .308<  0.0)1 . 00798 1 <-.00841  > 

. 006401  . 00793 ) ( 186 .9)  < .00848) 

.00143£<  0.0)  < . 001 432 > 

- .0014321 0. OX 0.0)  <-. 001438 > 

.001430<-  .0317)  <-.454E-4> 

.00841 < .373E-9)  <.906E-18) 

-.00240C  .26952.23  .64752.14)  <-.01800^ 

• 00£40<  .00600)1  .38251  .99  .643*1.888)  <.765E-4> 


MUMEP8T0RS  : 

.0006381 .1333)1  .504*33.3  16.76*88.7)  <.0971) 

-.0003961  0.0X-6.41  ) < .00854) 

-.0006561  159.1)1  .9665.964  .9305.850)  <-.0971) 

.0006581-  .000609X159.1)1  .991  * .94  3 .934*  .1266) 

■580E-10) 

-.000657(139.1)1  .966  5.964  .930  5.850)  <-.0970) 

.765E-4I - .888 '(39.3)  <-.000867 

.016140-. £68X39 .3)  <-.1830-' 

- .016141  .00799X-. 888X39. 3)  < .001461  > 

. J44E-4I -.864X48.5X460.)  <-.1776) 

.7 65E— 4<-.888X 39.3)  -.000867) 

-.76^E-4i  0.0X-. 868X39. 3)  < .000867) 

•356E-4<-. 277X56. 3)  <-.000666 

. 00018691 - .866  X 39 .3)  <-.00146t  > 

-.344E-4I  1 59.1  )t  .9665.964  .9305.830)  <-.00308) 

. 344E-4I  1 .937X159.1)1-.  11761. 371  - .04365.366) 

<.  .0014*9  ) 
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STABILITY  DERIVATIVES  AND  TRANSFER  FUNCTION  FACTORS 
(STOL  COMMUTER,  LANDING) 


C-EDMETRY: 


VT 

8mMMh 

119 .23 

-”7  .999 

8 

PHD 

Mfli;  h 

: : 10 

2J 

1 1 IF  .5 

. 002377 

. 1 0590 

45.20 

-2.043 

s £ 

c 

HEIGHT 

IV 

ALTITUDE 

420.0 

6.500 

11001. 

23470. 

0.0 

MOM-PIMEM:  IOH8L 

per  ivrtive: 

■S 

CL 

CLft 

C LmD 

CLM 

1 .*19  04 

5.251 

—9 . 190 

0.0 

CM8 

CMHl' 

CMG 

CMM 

-1 .4991 

-6 . 060 

-32.84 

0.0 

CD 

COR 

1 DM 

TM 

tdth 

.17212 

.4004 

0.0 

-8320. 

42 . 07 

CD  PE 

CLDE 

CM  PE 

0.0 

.4>-  70 

- 1 .8219 

D I MEMS  1 HMHL  DERIVATIVES  ! 

XU 

::u  ♦ 

XU 

TU 

- . 0594c 

- .07476 

.19  31 

-.02180 

zu 

71 J ♦ 

rwr 

2G 

- . 5 367 

- .52  ?£’ 

.0  881 

- .9361 

MU 

MU  ♦ 

MUD 

MU 

.0001 055 6 

0007546 

.002722 

- . 02429 

MHT' 

MR 

MG 

- . 'Jcl  1 9 

~c  .671 

-1 .7443 

:;iiE 

ZI'E 

MDE 

0.0 

— M _ V " C| 

- 3 .*2 0 

XDTH 

2DTH 

moth 

.06670 

-.0-731 

. 00  7*62 

PEMnMIhHTDK1.: 

C .1271  ; . 

11  .0?95! 

. 3 06  3C  . • 

0452.16  1 .92411  ,?37] 

* .453 

DE  riUMFBHTnE  : 

M<  U-PE  > -t  -99Ti(2.66)<-18.29>  <99.5 

M'M-DE’  -10.  4<43.6'C  .13298.372  . 044”  !.  3*9  3 < -92  .5 > 

fl<  THE-DF ' -3.44<  .222 '«  .763)  -.540> 

Nipp-DEi  I 0 . 34  ‘ . 02  01  > ' - 4 .33M8 . 07)  07.2* '> 

M<8r-PE  -1(1.  4<-4.44  ><8.15it  .349!  .0728  . 0234  ! . 06*  3 3 

1 . 99*  > 

M>  HP-PE  > 10.2-(  .049*. ''-4.36X8.10)  -17.61) 


PTH  Mi'MER^TieS* 

rt.U-DTM.  .<)>•►  • .t£9”o[  .65712.08  t . 366  ! ! .569  3 < . 0486  > 
Mi'-PTM  - r.l  1 ( .24,"!  . !22  .0994!.  307  3 -.0668> 


TAB  IE  B-4  (Continued) 


PE  PTH  COl'FLIMG  Ml  'MERWTnDT  : 

M<U-DE.U-BTH)  .896<-. £03X43.8)  <-7.96> 

M<  U-DE - THE-PTM)  .3031 .739)  < .££4 
H<U-PE  BD-DTH>  - .896<-4  .82  > <8  . : 1 ) <34.4> 
h<U-PE • HZ-PTH > .898*  .0818X-4. 72X8. 39)  <-.753> 

MUI-PE^HP-BTH)  -.8981-4. 68X8. 31)  <34. £> 
NIW-PE'THE-PTH’  -.8881.667)  <-.1879> 

M*W-PE  PP-PTH  * -33.3*  .887)  <-88.8> 

M*W-PE. HZ-PTH < 33.31 .0885)1 .667)  <.63£ 

N<W-PE  HP-PTH  ■ . 0937 < .787) 1-311.)  <-8£.9> 

M«  THE-PE  - DP-PTM)  -.8881.667)  1-.1879) 
NlTHE-PE'HZ-PTH)  .£8810.0)1.667)  <.1879> 

M( THE-PE  HP-PTH)  -.2491.657)  <-.1635> 
NIPP-PE'HZ-PTH ' - .9481 .887 ) <-.632> 

h*  PD-PE -HP-PTH)  - . 0937*  -4.88X8.31)  <3.60> 

Ml HZ- PE 'HP-PTH)  .09371-  .80“  ) <-3  .75X7  .85)  < .550' 


UG  NUHEB9T0RS » 
M<  IJ-UG  ' 
M'W-l.'G' 
h<  THE-UG) 
N<  PP-IJ6) 

MIHZ-IJG) 

HIHP-UG) 


.074818. 04)C  .881  51  .788  1.172(1.8613  < .453) 

.544*  -.883E-6X-  .00308)1 1 .577)  < .587E-10) 

-.00884*  0.0X8.84)  <-.01398) 

-.544* -.745E-8  t .59851 .740  1 .03051.4083 

- .575E-9) 

.544<  .671E-7X  .0889)1  .58851.787  1 .01551 .397  3 

<- .454E-1 0> 

-.549*  .0884)1  .59851  .748  1.03851.4033  <-.0473) 


WG  MUMERHTORS : 

M'  1.1- MG) 

M<  U-UG  ■ 

M<  THE -MG  > 
M>  PP-MG ■ 

M*  83 -MG ) 
fi.  HP-MG  . 


- . 1907*  - . 373E-8 >C  .548  51.371  .743  51.158  3 

<- . 181  IE-9 

-.  04041-5. 88X20. 1)1  .14775.309  .04565.3053 

.453) 

- .01  19r.  • o . ox  - . 0896  ) < - .7*0)  <-  .000801 

. 0404*  3.58*  -84  .8  •!  .1457  5 . 383  . 0589  5 . 359  ] 

-.453) 

.830<0. 0>«4. 19  1 .13155 .368  .04765.359]  .455) 

.0408*  3 .57)*  -23.7  )t  .13185.364  .04805.361  3 

- .450. 
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TABLE  B-k  (Continued) 

I'E  < 06  COWLING  MIJMES6T0R1  : 

M1U-BE  W-IJG  ■ ,773< 1 .940X  41  .7)  62.5? 

M<U-IlE  TME-UG?  .26112.26)  .590? 

Me.  U-FE  - M-U6  ‘ -.7731-1.730X5.42)  <7. 25) 

McU-I'E  83-06  > .773<  .227X-2. 05X5.51)  <-l  .986> 

M<  U-PE  HB-UG  > -.7691-1  .730X5.42)  <7.21> 

M < M-BE - THE -UG  > 1 . 924 t 0 . 0 > < t . 924 > 

MlW-BE- BB-IJ6)  227  . > - . 161 9E-9 ) <-.368E-7> 

MlW-BE  87-06)  -227. e .489E-9X .0285)  <-.316E-8> 

fCM-PE  HB-UG)  .0808<  .0284X2840.)  <8.53  : 

M<THF.-BE  BB-U6)  1 -924c  .622E-9 > <.H98E-8> 
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TABLE  B-1+  (Concluded) 
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TABLE  B-5.  STABILITY  DERIVATIVES  (JUMBO,  TAKEOFF) 
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TABLE  B-6.  STABILITY  DERIVATIVES  (LIGHT  TRANSPORT,  TAKEOFF) 
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TABLE  B-8.  STABILITY  DERIVATIVES  (STOL  COMMUTER,  TAKEOFF) 
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APPENDIX  C 


DERIVATION  OF  SIMPLIFIED  PILOT- VEHIC 125  DYNAMICS 


This  appendix  develops  a simplified  description  of  airspeed- flight  path 
dynamics  useful  for  analyzing  wind  shear  effects.  The  ultimate  objective 
is  to  directly  express  airspeed  and  flight  path  excursions  for  the  closed 
loop  pilot- vehicle  in  terms  of  essential  parameters,  i.e.,  aircraft  stability 
derivatives  and  pilot  crossover  frequencies  (loop  tightness).  The  general 
form  is  set  up  for  both  CTOL  and  STOL  piloting  techniques,  but  only  the  CTOL 
case  is  fully  developed.  The  STOL  case  would  follow  in  the  same  manner, 
however . 

The  approximate  pitch- constrained  equations  of  motion,  without  flight 
path  or  airspeed  regulation,  are: 


s - X. 


u xws 


V 

) 

Ah 

J 

_ 

L 

OCa  — g)  X5t  -s  — Xw 
— Za  ~ Z&T  0 Zw 


0 

&T 


u, 


g 


The  elements  of  the  pilot  model  are: 
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or 
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Considering  the  CTOL  case  as  an  example  and  substituting  the  above  pilot 
model  into  the  pitch  constrained  equations: 
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The  airspeed  and  altitude-due- to-gust  transfer  functions  are  easily 
obtained  by  substituting  appropriate  columns  in  the  matrices  and  solving 
for  determinants  (see  Ref.  25  for  details).  For  example,  in  the  case  of 
altitude,  pilot  gains  are  set  according  to  an  airspeed  loop  tightness 
(crossover  frequency,  a>Cu)  and  flight  path  loop  tightness  (crossover  fre- 
quency, . A good  first  order  approximation  is 


and 


-KdZa  = “cd 
Ku^&t  = ^u 


Therefore,  the  equations  of  motion  are  approximately; 
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The  characteristic  determinant  of  the  above  equations  of  motion  will  factor 
into  a dominant  airspeed  mode,  l/Tu,  and  an  oscillatory  dominant  flight  path 
mode  characterized  by  a damping  ratio  £d  and  natural  frequency  o)d.  A good 
approximation  for  the  speed  mode  of  CTOL  aircraft  appears  to  be: 


Tu 

This  permits  an  approximate  solution  of  the  flight  path  mode. 
To  summarize, 
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